{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 4. Noise Reduction in NDVI Timeseries\n", "*Written by Men Vuthy, 2021*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Import packages**" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Read data\n", "DF_NDVI = pd.read_csv('output/2/timeseries/timeseries_ndvi.csv')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Remove outlier using Hamper Filter**" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "def hampel_filter_forloop(input_series, window_size, n_sigmas=1):\n", " \n", " n = len(input_series)\n", " new_series = input_series.copy()\n", " k = 1.4826 # scale factor for Gaussian distribution\n", " \n", " indices = []\n", " \n", " # possibly use np.nanmedian \n", " for i in range((window_size),(n - window_size)):\n", " x0 = np.median(input_series[(i - window_size):(i + window_size)])\n", " S0 = k * np.median(np.abs(input_series[(i - window_size):(i + window_size)] - x0))\n", " if (np.abs(input_series[i] - x0) > n_sigmas * S0):\n", " new_series[i] = x0\n", " indices.append(i)\n", " \n", " return new_series, indices" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Prepare input for hamper filter\n", "HF_input = np.array(DF_NDVI.T)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "result = []\n", "\n", "for i in range(len(DF_NDVI.T)):\n", " res = hampel_filter_forloop(HF_input[i], window_size=5, n_sigmas=2)\n", " result.append(res[0])" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "NDVI_HF = pd.DataFrame(result).T" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "NDVI_HF.columns = DF_NDVI.columns" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Smoothen data using Moving Average method**" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "NDVI_MA = []\n", "\n", "for i in range(len(NDVI_HF.columns)):\n", " column = NDVI_HF.iloc[:,i:(i+1)]\n", " mvg_avg = column.rolling(window=5).mean()\n", " NDVI_MA.append(mvg_avg)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "NDVI_MA = np.array(NDVI_MA)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(1, 69, 192289)" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# In order to export, write requires an array of shape (band, row, col). so we have to reshape the axis\n", "NDVI_MA = np.moveaxis(NDVI_MA, [0, 1, 2], [2, 1, 0])\n", "NDVI_MA.shape" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "NDVI_MA = pd.DataFrame(NDVI_MA[0])\n", "NDVI_MA.columns = NDVI_HF.columns" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Plot result**" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "i = 5000 # column number or pixel numer\n", "\n", "# Create data frame\n", "j = i+1\n", "noised_data = NDVI_HF.iloc[:,i:j]\n", "smooth_data = NDVI_MA.iloc[:,i:j]\n", "df = pd.concat([noised_data, smooth_data], axis=1)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'ndvi_2010-2012')" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAskAAAE/CAYAAAC0Fl50AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAACnDElEQVR4nOzdd3iUVfbA8e+dmUx6740kQCDUoFTpiCDYe2+oq66r6+66Rbe5v3V3davrWtfe64KiWBEpAlJCbwFCCgnpvWcyM/f3xzsTQkiZlkbu53l8gJn3fecGSebMec89R0gpURRFURRFURTlJF1/L0BRFEVRFEVRBhoVJCuKoiiKoihKBypIVhRFURRFUZQOVJCsKIqiKIqiKB2oIFlRFEVRFEVROlBBsqIoiqIoiqJ0oIJkRVEURVEURelABcmKoiguEEJIIcRIB477tRDipb5Yk6IoiuI5KkhWFEXpRVLKv0gp7+zuGCHEKCHESiFEmRCiUgjxlRBidIdjfiqEKBZC1AghXhFCeLd77j4hRIYQokUI8Von118ohMgUQjQKIdYKIZK6WcsMIcRq2zrKhBAfCiFi2z0vhBB/FUJU2P77mxBCtHv+USHEPiGEWQjxhw7XvlAIsVEIUW37Wl4UQgR293ejKIrSX1SQrCiK0v9CgE+A0UA0sA1YaX9SCHE+8BCwEEgGhgP/1+78QuBPwCsdLyyEiABWAL8DwoAM4P1u1hIKvGB7nSSgDni13fN3AZcB6cBE4CLg7nbPZwG/BD7r5NrBtnXGAWOABODv3axFURSl36ggWVEUxUYIkSuE+LkQYq8tY/u+EMLH9twvhBBFQohCIcTt7c6ZYcuK6ts9drkQYq/t938QQrzV3etKKbdJKV+WUlZKKVuBJ4DRQohw2yG3Ai9LKQ9IKauAR4Hb2p2/Qkr5MVDRyeWvAA5IKT+UUjYDfwDShRBpXazlC9uxtVLKRuBpYFa7Q24F/imlLJBSngD+2WEtr0spv0ALrjte+x0p5ZdSykbb1/Fih2sriqIMGCpIVhRFOdU1wBIgBS1TepsQYgnwc2ARkAqcZz9YSrkFaADObXeNG4B33FjDXKBYSmkPescBe9o9vweIbhdEd+eUc6WUDcAx2+OOruVAV9ez/d7Ra/V0bUVRlAFDBcmKoiin+o+UslBKWQl8CkxCC5xflVLutwWZf+hwzrvA9QC2GtsLbI85TQiRADwD/KzdwwFATbs/23/vSD1vx3Pt5/d4rhBiIvB74Bc9rCWgfV2yI4QQi9Cy0r935jxFUZS+ooJkRVGUUxW3+30jWlAYB+S3ezyvwznvAFfYNtNdAeyUUnY8pkdCiEjga+BZKWX7ILseCGr3Z/vvTytp6ETHc+3n1wkhhgkh6u3/dVjLSOAL4AEp5Xc9rKVeSikdWIv92jPQ/s6uklIecfQ8RVGUvqSCZEVRlJ4VAYnt/jys/ZNSyoNogfNSXCy1EEKEogXIn0gp/9zh6QNoG+Xs0oGSduUY3TnlXCGEPzACrU75uJQywP5fu2OSgG+AR6WUbzqwFodLJoQQZ6FtUrxdSrnG0fMURVH6mgqSFUVRevYBWm3yWCGEH/BIJ8e8A/wYrc72Q2cuLoQIAr4CNkkpH+rkkDeAO2yvHwr8Fnit3fkG2wZDPaAXQvgIIQy2pz8CxgshrrQd83tgr5Qys4u1xAPfAs9IKZ/vYi0/E0LECyHigAc7rMXL9jo6wGBbi9723HjgS+B+KeWnjv3tKIqi9A8VJCuKovTA1q3h32jBY5bt147eBeYD30opy518icuBqcCy9uUPQohhttf/EvgbsBYtY53HqYH6b4EmtDZxN9l+/1vbuWXAlcCfgSpgOnBdN2u5E63F3CNdlGL8F61Wex+wH63V23/bPf+i7fWvB35j+/3NtuceBCKBl9tdW23cUxRlQBJOlJEpiqIoiqIoypCgMsmKoiiKoiiK0oEKkhVFUfqAEOLGDqUUqtxAURRlAFPlFoqiKIqiKIrSgcokK4qiKIqiKEoHhp4P6XsREREyOTm5v5ehKIqiKIqinMF27NhRLqWM7Oy5ARkkJycnk5GR0d/LUBRFURRFUc5gQogup6OqcgtFURRFURRF6UAFyYqiKIqiKIrSgQqSFUVRFEVRFKWDAVmTrCiKoiiKorintbWVgoICmpub+3sp/c7Hx4eEhAS8vLwcPkcFyYqiKIqiKGeggoICAgMDSU5ORgjR38vpN1JKKioqKCgoICUlxeHzHCq3EEIsEUIcFkJkCSEe6uT5YCHEp0KIPUKIA0KIZe2eyxVC7BNC7BZCqJYViqIoiqIofaC5uZnw8PAhHSADCCEIDw93OqPeYyZZCKEHngEWAQXAdiHEJ1LKg+0O+xFwUEp5sRAiEjgshHhbSmmyPb9ASlnu1MoURVEURVEUtwz1ANnOlb8HRzLJ04AsKWW2Leh9D7i0wzESCBTaCgKASsDs9GoURVEURVGUM0p1dTVXXXUVaWlpjBkzhu+//55f/OIXpKWlMXHiRC6//HKqq6sBMJlMLFu2jAkTJpCens66devarvPuu+8yYcIEJk6cyJIlSygv1/KveXl5LFy4kIkTJzJ//nwKCgo8sm5HguR4IL/dnwtsj7X3NDAGKAT2AQ9IKa225yTwtRBihxDiLjfXqyiKoiiKogwiDzzwAEuWLCEzM5M9e/YwZswYFi1axP79+9m7dy+jRo3iscceA+DFF18EYN++faxevZoHH3wQq9WK2WzmgQceYO3atezdu5eJEyfy9NNPA/Dzn/+cW265hb179/L73/+ehx9+2CPrdiRI7iw/LTv8+XxgNxAHTAKeFkIE2Z6bJaU8G1gK/EgIMbfTFxHiLiFEhhAio6yszJG1K4qinPFaLVbWHS7t72UoiqK4pLa2lg0bNnDHHXcAYDQaCQkJYfHixRgMWtXvjBkz2rK/Bw8eZOHChQBERUUREhJCRkYGUkqklDQ0NCClpLa2lri4uNPOWbBgAStXrvTI2h0JkguAxHZ/TkDLGLe3DFghNVlADpAGIKUstP1aCnyEVr5xGinlC1LKKVLKKZGRnY7QVhRFGXI+yMjntle3s6+gpr+XoiiK4rTs7GwiIyNZtmwZZ511FnfeeScNDQ2nHPPKK6+wdOlSANLT01m5ciVms5mcnBx27NhBfn4+Xl5ePPfcc0yYMIG4uDgOHjzYFninp6ezfPlyAD766CPq6uqoqKhwe+2OtIDbDqQKIVKAE8B1wA0djjkOLAS+E0JEA6OBbCGEP6CTUtbZfr8Y+KPbq1YURRki1hzSssgHi2qYkBDcz6tRFGWw+r9PD3CwsNaj1xwbF8QjF4/r9hiz2czOnTt56qmnmD59Og888ACPP/44jz76KAB//vOfMRgM3HjjjQDcfvvtHDp0iClTppCUlMTMmTMxGAy0trby3HPPsWvXLoYPH87999/PY489xm9/+1v+8Y9/cN999/Haa68xd+5c4uPj27LU7ujxClJKsxDiPuArQA+8IqU8IIS4x/b888CjwGtCiH1o5Rm/klKWCyGGAx/ZdhQagHeklF+6vWpFUZQhoMlkYVOWtjHlUFFdP69GURTFeQkJCSQkJDB9+nQArrrqKh5//HEAXn/9dVatWsWaNWvauk8YDAaeeOKJtvNnzpxJamoqu3fvBmDEiBEAXHPNNW3XiYuLY8WKFQDU19ezfPlygoPdTyo4FGZLKT8HPu/w2PPtfl+IliXueF42kO7mGhVFUYakzcfKaTFb8TboyCz2bAZIUZShpaeMb2+JiYkhMTGRw4cPM3r0aNasWcPYsWP58ssv+etf/8r69evx8/NrO76xsREpJf7+/qxevRqDwcDYsWMpLCzk4MGDlJWVERkZyerVqxkzZgwA5eXlhIWFodPpeOyxx7j99ts9snY1cU9RFGWA+uZQKf5GPUvGx/JtZglSStXzVFGUQeepp57ixhtvxGQyMXz4cF599VWmTp1KS0sLixYtArTNe88//zylpaWcf/756HQ64uPjefPNNwEtW/zII48wd+5cvLy8SEpK4rXXXgNg3bp1PPzwwwghmDt3Ls8884xH1i2k7Nioov9NmTJFZmSo4XyKogxdUkpmPLaGsxJDOWdEOI98coCtv15IdJBPfy9NUZRB4tChQ23ZVqXzvw8hxA4p5ZTOjndoLLWiKIrStw4U1lJS28K5Y6JIiwkE4FCRKrlQFEXpKypIVhRFGYC+zSxFCFgwOoq0GK3tfGax2rynKIrSV1RNsqIoygC0JrOU9IQQIgO9AYgL9uGwCpIVRVH6jMokK4qiDDBldS3sya9mYVpU22OjYwJVuYWiKEofUkGyoijKALM2Uxsgcu6Yk0FyWmwQx8rqMZmt/bUsRVGUIUUFyYqiKAPMmswSYoN9GBsb1PZYWkwgrRZJdnl9P65MURRl6FBBsqIoygDSYrbw3dFyzk2LOqUn8hhbwKzqkhVFGWyqq6u56qqrSEtLY8yYMXz//ff84he/IC0tjYkTJ3L55ZdTXV0NgMlkYtmyZUyYMIH09HTWrVvXdp13332XCRMmMHHiRJYsWUJ5uTaRNC8vj4ULFzJx4kTmz59PQUGBR9atgmRFUZQBZGt2JY0mCwvblVoApET446UXajy1oiiDzgMPPMCSJUvIzMxkz549jBkzhkWLFrF//3727t3LqFGjeOyxxwB48cUXAdi3bx+rV6/mwQcfxGq1YjabeeCBB1i7di179+5l4sSJPP300wD8/Oc/55ZbbmHv3r38/ve/5+GHH/bIulWQrCiKMoB8m1mKj5eOmSMiTnncS69jZFSgGk+tKMqgUltby4YNG7jjjjsAMBqNhISEsHjxYgwGrcnajBkz2rK/Bw8eZOHChQBERUUREhJCRkYGUkqklDQ0NCClpLa2lri4uNPOWbBgAStXrvTI2lWQrCiKMkBIKVmTWcKsERH4eOlPe35MTCCZKpOsKMogkp2dTWRkJMuWLeOss87izjvvpKGh4ZRjXnnlFZYuXQpAeno6K1euxGw2k5OTw44dO8jPz8fLy4vnnnuOCRMmEBcXx8GDB9sC7/T0dJYvXw7ARx99RF1dHRUVFW6vXfVJVhRFGSCOltaTX9nEPfNGdPp8WmwgK3adoLrRRIifsY9XpyjKoPbFQ1C8z7PXjJkASx/v9hCz2czOnTt56qmnmD59Og888ACPP/44jz76KAB//vOfMRgM3HjjjQDcfvvtHDp0iClTppCUlMTMmTMxGAy0trby3HPPsWvXLoYPH87999/PY489xm9/+1v+8Y9/cN999/Haa68xd+5c4uPj27LU7lBBsqIoygCx5pCt9VtaVKfPj243eW/G8PA+W5eiKIqrEhISSEhIYPr06QBcddVVPP64Fli//vrrrFq1ijVr1rRtVDYYDDzxxBNt58+cOZPU1FR2794NwIgRWhLhmmuuabtOXFwcK1asAKC+vp7ly5cTHBzs9tpVkKwoijJAfJtZwri4IGKDfTt9fkxMIACZRbUqSFYUxTk9ZHx7S0xMDImJiRw+fJjRo0ezZs0axo4dy5dffslf//pX1q9fj5+fX9vxjY2NSCnx9/dn9erVGAwGxo4dS2FhIQcPHqSsrIzIyEhWr17NmDFjACgvLycsLAydTsdjjz3G7bff7pG1qyBZURRlAKhqMLEjr4r7Fozs8pjIQG/C/I1kqjZwijKoVDaYqG1qJTnCv7+X0i+eeuopbrzxRkwmE8OHD+fVV19l6tSptLS0sGjRIkDbvPf8889TWlrK+eefj06nIz4+njfffBPQssWPPPIIc+fOxcvLi6SkJF577TUA1q1bx8MPP4wQgrlz5/LMM894ZN0qSFYURRkA1h8pwyrh3DHRXR4jhCAtJlAFyYoyyDz+xSG25lSy/hcL+nsp/WLSpElkZGSc8lhWVlanxyYnJ3P48OFOn7vnnnu45557Tnv8qquu4qqrrnJ/oR2o7hYDVGF1E40mc38vQ1GUPrIms5SIACMT47uvoxsdE8jh4jqsVtlHK1MUxV255Y0cr2ykxWzp76UoTlBB8gCUX9nIef9azxOrj/T3UhRF6QOtFivrDpeyYHQUOp3o9tgxMUE0tVo4XtnYR6tTFMVdhTVNSAlF1c39vRTFCSpIHmCklPz6o300mizsLajp7+UoitIHMnKrqGs2nzZlrzNpsbbNe2qoiKIMClarpKRWC44Lqpr6eTWKM1SQPMAs33mC746WExFg5EhJHVKqW6qKcqb7NrMEo17H7NTIHo9NjQpEJ1B1yYoySJTXt9Bq0d7LC6r6/g6QiiM0rvw9qCB5ACmra+HRVQeZnBTKD+ePpKqxlbL6lv5elqIovWxNZinTh4cR4N3zXmpfo57kcH81eU9RBonCmpMlFn2dSfbx8aGiomLIB8pSSioqKvDx8XHqPNXdYgD5w6cHaDJZ+OuVEyit1YLjI8X1RAU69z9VUZTBI6e8geyyBm6ZkeTwOWmxgRwsVOUWijIYFFVrgbEQkN/HmeSEhAQKCgooKyvr09cdiHx8fEhISHDqHIeCZCHEEuBJQA+8JKV8vMPzwcBbwDDbNf8hpXzVkXMVzeqDJXy2t4gHF41iZFRg28jZIyV1zE6N6OfVKYrSW77NtE/Z67r1W0dpMUF8sb+YRpMZP6PKdSjKQGbPJI+ODuzzTLKXlxcpKSl9+ppnkh7LLYQQeuAZYCkwFrheCDG2w2E/Ag5KKdOB+cA/hRBGB88d8mqbW/ntx/tIiwnk7nnauMWIAG/C/bW6ZEVRzlxrDpWQGhXAsHC/ng+2SYsJREo4UlLfiytTFMUTiqqb8PXSMz4+uF9qkhXXOVKTPA3IklJmSylNwHvApR2OkUCg0AZvBwCVgNnBc4e8x7/IpKyuhb9eORGj4eT/ktToAA6rIFlRzli1za1sy6lkYTcDRDqTFhMEaOOpFUUZ2AprmogN8SEx1I+S2hbVK3kQcSRIjgfy2/25wPZYe08DY4BCYB/wgJTS6uC5Q9qW7Are2Xqc22elkJ4Ycspzo6MDOVKsOlwoyplqc1YFZqtkweieu1q0lxDqi79RrzpcKMogUFjdTFywLwmhvm1/VgYHR4Lkzjrbd4zazgd2A3HAJOBpIUSQg+dqLyLEXUKIDCFExlApMG9utfDwin0khvnys8WjTnt+VEwgDSYLJ6pVX0VFORNtya7Ax0vHWcNCnTpPpxOMjgnkkMokK8qAV1TTRGywT1uQnK8GAQ0ajgTJBUBiuz8noGWM21sGrJCaLCAHSHPwXACklC9IKadIKadERjqXVRmsnlxzlJzyBh67fGKnm29GR2tDA46qukNFOSNtya5gclLoKWVWjkqLDSJT3WlSlAGt1WKltK6F2BBfEsK0fQdqoMjg4chP5u1AqhAiRQhhBK4DPulwzHFgIYAQIhoYDWQ7eO6QtP9EDS9syOaaKQlddq9ItQXJqi5ZUc481Y0mDpfUMSMl3KXz02ICqWlqpaRW9VJXlIGqpLYZKSEu2IeYIB8MOqE27w0iPfYOklKahRD3AV+htXF7RUp5QAhxj+3554FHgdeEEPvQSix+JaUsB+js3N75UgYPs8XKr5bvJczfyG8u6LrZR7CvF7HBPhxRdYeKcsbZmlOJlDB9uKtBsrZ571BxLTHBqpe6ogxERbb2b7Ehvuh1grgQX5VJHkQcarAppfwc+LzDY8+3+30hsNjRc4e6F7/L4UBhLc/deDbBfl7dHpsaHagyyYpyBtqaXYm3QUd6YrBL54+O0e40ZRbVsWB0lCeXpiiKhxTa9hTF2T7IJoT6qkzyIKLGUvcxKSUvbDjGgtGRLJ0Q2+Pxo6MDyCqtx2JVdYeKcibZkl3B2cNC8TboXTo/2NeL+BBfMovV5j1FGajaZ5JBC5LzVSZ50FBBch87XtlIVWMr5411rC/qqOhAWsxWjqvdsIpyxqhpbOVQcS0zXCy1sBsdE8hhVY6lKANWUXUTgT4GAry1G/cJoX6U1bXQ3Kp6JQ8GKkjuY3sKagBITwhx6Hj7LVX1RqgoZ45tuVo98ozhYW5dJy0mkKzSekxmq4dWpiiKJxXWaD2S7RLDtN+r1q6DgwqS+9je/GqMBl1b8NuTkVEBCIEaT60oZ5Ct2RUYDbrTBgg5Ky02CLNVcqxMtYlUlIGoyDZtzy4hVLWBG0xUkNzH9p6oYWxsEF56x/7q/YwGEkP91OY9RTmDbMmp4OxhIfh4uVaPbDfGvnlP1SUryoBUVN1MbLtMshooMrioILkPWayS/SdqSE9wbjf7KNt4akVRBr+aplYOFtYy3cX+yO0lR/hj1OvUeGpFGYCaWy1UNJjaOlsARAX64KUXKpM8SKgguQ8dK6un0WRhooP1yHajYwLIKW9QdYeKcgbIyK3EKnF70x6Al17HyKgAMotUkKwoA01xh84WAHqdID5EtYEbLFSQ3If25FcDON0XdVR0IGarJKe8oRdWpShKX9qaU4lRr+OsYSEeuV5abKAqt1CUAaiw5tQeyXYJoX4qkzxIqCC5D+0tqCHA28DwiACnzmvrcKHqkhVl0NuSXcEkD9Qj242JCaKktoWqBpNHrqco/ammsZWnvz2K2TL475wWVp+eSQb7QBEVJA8GKkjuQ3sLqhkfH4ROJ5w6LyXCH71OqLpkRRnk6ppb2X+ihhkp7rV+a69t8p76+aCcAT7bV8Q/vj7CnoLq/l6K24psbd5iT8sk+1Je30KTSfVKHuhUkNxHTGYrh4rqHO6P3J63QU9KhL/KJCvKIJeRW+WxemS7tNhAxokcYr66C2oLPXZdRekPeZVaWWFO+eCv2S2saSbc33jaXSN7G7gT1YP/azzTqSC5j2QW12KyWJ3etGc3OjpQ9UpWlEFuS04FXnrBWcNCPXbNSD89T3i/QErpN/Dm5dBQ4bFrK0pfs7dGyz0D9uB07JFsZx8oosZTD3wqSO4j9kl7E51s/2Y3KjqQ45WN6vaMogxiW7IrmZQYgq/RM/XIACLjFUaRx0e+l0NlDrx9JTSrjXzK4JRXoQXJORVnQJDcoUeynRooMnioILmP7M2vJszf2NZI3FmjYwKQErJK1WQtRRmM6lvM7D9R45H+yCcvWgrf/oljgdP4df21WK56DYr2wns3QKt6A1YGFyklx+1BctngD5ILa5pO62wBEBngjVGvU23gBgEVJPeRvQU1TEwIRgjnNu3ZjYpWHS4UZTDLyK3EYpUerUdm9SPQ2sThs39LU6uVE1Hz4fL/Qu5G+HAZWFo991qK0suqGlupazFjNOjIrWhASunytaSUtJj7785rfYuZumbzaZ0tAHQ6QXyoLwWV6oPsQKeC5D7QaDJztLTO5XpkgKRwf4wGnapLVpRBamtOJQad4OykEM9c8PgW2PMOzLyPwISxgFYDycSr4cJ/wJEv4ON7wTr4W2kpQ0OercRiekoYjSYLZXUtLl9r5e5CpvzpG7eu4Y6uOlvYaW3gVCZ5oFNBch/Yf6IWq8TpcdTt6XWCkZEBHFZtnhRlUNqSXUF6Ygh+RoP7F7OY4bOfQ1A8zP0FMUHaG3FxrdaXlal3wrm/g30fwBe/ADcycorSV47bNu3NGxUJQLYbm/e251ZS12xm+c4Cj6zNWYW2aXtxnWSSQQ0UGSxUkNwH9tr6PbqTSQatH+qRkjrtDa+xEsqzoKna7fUpitK7GlrM7CuoYbqn+iNnvAIl++D8P4PRn2hbtqrEHiQDzHkQZt4P21+Cb//kmddVlF5kr0e2B8nudLiw7995f3u+W2UbrnIkk1zRYKLRZO7LZSlO8kBKQ+nJnoIa4oJ9iAz0dvyknW9om3IayqGxHBrK+E3JCSzNZchH6xHWdt9YQfEQNcb231jt14jRYPTz/BejKIrTduRVYfZUPXJ9mRb0Dp8PYy8DINDbgJ9RT3FNu1vLQsCiR6G5Br77B/gEw6wfu//6itJL8iobiQ7yZnhkAEa9zq0OF1ml9YT5G8kpb2BrTqVn9wI4oLCmGSEgOqjrIBngRFUTqbY9R8rAo4LkPrC3oNr5LPLq30NTFRgDwT8c/CORIcNYWxPPvPSxxMYmgl8Y1BVB6SEoPQg534HF/iYpICxFC5oXPgKRozz9ZSmK4qCtORXodYLJSR7oj/zNI9DaCEv/rgXCgBCCmCCfUzPJ2hNw0b+hpQ5W/w68A2HKMvfXoCi94HhFI0lh2oTZxDBflztcVNS3UNFg4sFFo3jhu2ze23a8z4PkouomogK98dJ3fsPe3gYuv6pRBckDmAqSe1l1o4m8ikaunZro3Ik//B58Q8DrZD1Tc2UjD/1tLX+Jn8AN04edfo7FDFU5WsBsD5yPrYVVP4HbPmt7Q1UUpW9tya5kYkIw/t5u/sg9vhV2vw2zfnLaB9+YYJ+TNcnt6fRw+QtaoLzqJ7DrTZh2N4y7DAxO3N1SlF6WV9nAnFSt1CIlIoBcFzPJ9lKLCQnBXDYpnvcz8vm/xlaC/bw8ttaeFNV03iPZzj5QRNUlD2yqJrmX7bUNEXF6HHVQ7CkBMkB8iC/+Rn3XHS70BohIhbGXwvyH4Jo3tM07eZsgZ4MLq1cUxV2NJjN7C6rd749stcDnD7Zt1usoJsiH4ppOgmQAgxGufRuW/k0bNPLRXfDEOK1sQ42yVgaA5lYLJbUtDAvTMqwpEX7kVTRitTpfT5xVpgXJqdGBXDctEZPZyke7+nYDX2FNE3GdTNuziwzwxtugU0HyAOdQkCyEWCKEOCyEyBJCPNTJ878QQuy2/bdfCGERQoTZnssVQuyzPZfh6S+gLzWazDS3Otd30b5pb3y8650t7HQ6QWp0oHMdLs6+RXtTXfsXtcNdUfrBzrxqWi2SGcPd3LSX8QoU74PFfwLvgNOejg7Wyi26DCq8fGD63fCjbXDTCoifAhv+AU+Mhw9uhbzN6meE0m/snS2SwrUgOTnCnxazlaLO7o704GhJPX5GPXHBPoyLC2ZCfDDv9eEGPilll9P27ISw9UpWbeAGtB6DZCGEHngGWAqMBa4XQoxtf4yU8u9SyklSyknAw8B6KWVlu0MW2J6f4rml9y2rVXLNf7/nB284F+fvKahheIQ/wb6euc0zOjqQo6VOBMlePjDnZ5C/BY5965E1KIriOHs98pRkN4Lk+jL49lFImQfjLu/0kJggH8xWSUWDqftr6XQwciHc8B78eBeccy9kr4NXl8Lzc2Df/1xfp6K4yN7Z4mQm2R9wbfJeVmk9I6MC2oZ3XTctkcziOvbY7uz2turGVppaLV12trBTbeAGPkcyydOALClltpTSBLwHXNrN8dcD73picQPJ2sOl7D9Ry3dHy9l/wvFvNG3TnvtZZLtRMYGU15sor3eiQfpZN0NwIqx7TGWKFKWPbc2uZHx8MAHu1CN/8wcwNcAFf+9yb4F9F/1pm/e6E5aiZaZ/dggu/g9IKyy/Q9sErCh9KK8tk6wFx21Bsgt1yUdL6xgZdfJuyyXpcfh66Xlv23EPrLRnhTVa4NtVj2S7hFBf8itVJnkgcyRIjgfy2/25wPbYaYQQfsASYHm7hyXwtRBihxDiLlcX2p+klDy77hhxwT74G/W8vDHHofNKapspqW1xuz9ye6Ntu2Cdmrxn8NZ6phZsh6xvPLYWRVG619xqYXd+tXulFoW7YfdbcM6PIHJ0l4fF2LJWXdYld8foB5NvhR+sgeBh8OVD2kZgRekjxysaCPQ2EGrbXBcd6IOPl87pXsm1za2U1LaQGnWyY0SgjxcXTYzlkz2F1Lf0/r/romrte7CnTHJiqB9Vja19sibFNY4EyZ2lLbpKR14MbOpQajFLSnk2WrnGj4QQczt9ESHuEkJkCCEyysrKHFhW39mWU8mOvCrumT+Ca6cO49M9hQ69Ee3JrwYgPdGDmeRo7dPxEWcn7026EUKGqdpkRelDO49XYbJYmeHOpr19H4LeqH3Q7cZpU/dc4eULix+Fkv2w8zXXr6MoTsqrbCQxzK+tREKnEySH+5PjZJBs72zRPpMMcN20YTSaLKza0/sbVYucyCSD1itZGZgcCZILgPb9yxKArv6VXUeHUgspZaHt11LgI7TyjdNIKV+QUk6RUk6JjIx0YFl959l1xwj3N3L15ESWzUrGKiWvbc7t8bx9J2rQ6wRjYz0XJEcGehPi58XhknrnTjQYtR3xhTvhyFceW4+iKF3bkl2JTsCUZBf7I0sJmau0WmSf7n+ORAQY0Qknyy06M/ZSSJ6jdb5orOz5eEXxgOMVjW2b9uxSIvydziRn2d4bUzsEyWcPC2FUdADvbc/v7DSPKqxpxqATRAR032LRHiSrzXsDlyNB8nYgVQiRIoQwogXCn3Q8SAgRDMwDVrZ7zF8IEWj/PbAY2O+JhfeVA4U1rD9Sxu2zU/A16kkM82PJ+Bje2ZpHQw+3SPYU1DAqOhBfo95j6xFCMCo60LlyC7v06yE0GdapbLKi9IWt2RWMjw8m0MfFjbulh6AqF9Iu7PFQg15HVGA3beAcJQQs/as2qW/tX9y7lqI4wGKV5Fc1MqxDkJwc4c/xykbMFqvD18oqq8do0JEYduq1hBBcO3UYu/OrySyu9ci6u1JU3UR0kA96XfezCewDRdTmvYGrxyBZSmkG7gO+Ag4BH0gpDwgh7hFC3NPu0MuBr6WU7T/2RQMbhRB7gG3AZ1LKLz23/N733LpjBHgbuGlGUttjd8weTm2zmf/t6LrvopRS27TngdZvHY22BclOt7PRe8G8X0HRHjj8ucfXpSjKSc2tFnblV7s36SvzM+3X0UsdOjy6q4EizooeB1PugIyXoeSA+9dTlG4U1zbTapEkhfmf8nhKhD9mq3QqiDxaUseIyIBOA9QrzorHqNfx3rbezSYX1jQT30OpBWh3f3y8dGrz3gDmUJ9kKeXnUspRUsoRUso/2x57Xkr5fLtjXpNSXtfhvGwpZbrtv3H2cweL3PIGPt9XxI0zhp3Swm1yUihnDwvhlU05WLroSZpf2UR1YysTPViPbDcqJpC6ZnOPb4ZWq+STPYWU1rU7bsI1EDYC1j4GVsc/nSuK4pxNWeWYzFb3Nu1lroKEqRAY49DhMUHe7meS7Rb8Wivx+OJX6s6T0qvybB0sOiu3AOc6XBy1tX/rTKi/kfPHx7BiZ4HTMw+cUVTTRGw3g0TshBCqDdwApybudeO/G45h0Ou4Y3bKac/dOWc4eRWNfHOopNNz99iGiDg9ac8Bo2w/ALobKtJksnDfuzv58bu7eL/9p2a9Qcsml+zT3oAVRekVr23OJTrIu23MrtNqCqBoN4y+wOFTYoI8lEkG8AuDc38Lud/BwZU9H68oLurYI9ku2dYOztG65EaTmYKqptPqkdu7fmoitc1mvtxf7OJqu2e1Sop7GEndXkKoLwXVKpM8UKkguQsltc0s33GCqycnEBV4+ifCxWOjSQj15aXvsjs9f29BNUaDjtExgZ0+745RPbSBK6lt5toXvueL/cUYdILSug49lSdcBeGpWt9klU1WFI87UlLHd0fLueWcZLz0Lv6YPfyF9mvaRQ6fEh3sQ12zmUaTh1pKTV4G0RPg69+CSb2RK70jr7IRg06c1jItIsBIoLfB4Q4X2bbBI90FyTOGh5MU7se7vdQzubyhhVaL7HYkdXsJob4qkzyAqSC5Cy9vzMFstXLX3OGdPm/Q61g2K4XtuVXstrV6a29PQQ1jY4Ncf4PsRqi/kahAbw4Xn97hYl9BDZc8vZFjpfW8dMsUksL9Th88otPD/Ieg9CAcUhkiRfG0Vzfl4m3Qcf20Ya5fJPMz7cNs5CiHT2lrA+epkgudXtvEV5MPm//jmWsqSgfHKxpJCPXF0OH9UghBcoTjbeDs02hTo7sOknU6wTVTEtmaU0l2mZNdohxwskdyh0zyoVWw+hEwnzoRMyHUj+rGVuqaWz2+FsV9KkjuRE1jK29vyeOiiXFt0386c82UBAK9DacNF7FYJftP1JDuwUl7HY2OOX089Rf7irj6v5sx6HT874czWTgmmogAbyrqOxlTO+5yiEyDdY+DtfdqsxRlqKlqMLFiZwFXnB1PmL/RtYs0VWtlDmmOl1qAh3old5Q8S/t5sfHfUN377bOUoSevsoFhXbzXJkf4k+tgTfLRknoMOtHt+zbA1ZMT0OsE72d4/t+zvUfyKVnx5lr45D7Y9G94+yrtzzaJqsPFgKaC5E688X0uDSYLP5w/otvjAn28uH76MD7fV8SJ6pP/wI+V1dNosnh00l5H9jZwVqtESsnT3x7lh2/vZGxsEB//aBZjYoMAiAj07nyEtU6v1SaXZcKBj3ptnYoy1Lyz7TgtZivLZp2+l8FhWd+A1exUqQWcnLrndq/kjhY9qv26+neeva6iYOuR3KEe2S4lwp8TVU20mHtO5hwtrSc5wr/HO7hRQT4sTIti+Y4CWp1oL+eIQlsm+ZRBIluehaYqmP1TyNsEr14AdVpN9MleySpIHohUkNxBk8nCq5tzWTA6si3Q7M6tM5MBeG3TyWxyb0za62hUdADNrVayyur56fu7+cfXR7hsUhzv/GAGkYEnG5hHBnhT1lmQDDD2Mogaq7LJiuIhrRYrb36fx5zUiLa9Ay7JXAX+URA/xanTTo6m7uJ73lUhidob/IGPIOc7z15bGdKqG03UNptP62xhlxLhh1XiUJu0Y6X13dYjt3fdtETK602s6WLzvauKaprwNujaxmvTWAmbn9Y+8J73B7jhfajMhpcWQdkRNVBkgFNBcgfvbT9OZYOJexeMdOj4+BBfLpgQy3vb8ttqivYW1BDgbWB4hGPfrK6wvwHf+NJWPt5dyM8Xj+KJayfh43Xq4JJwfyN1zebO293odFptcsVROPhxr61VUYaKL/YXU1zbzLJZya5fxNwCR1drvZF1zv2I9jMaCPQxUFzTC1mpWT+G4GHw5UNg8dDGQOWMYjI7n5XN66KzhZ29w0VOefdBZIvZQm5Fg8NB8rxRUcQE+fCuh3smF9Y0Exfi2zZem03/BlO91ikGYOR5sOwzMDfBK4sJq9yFr5deZZIHKBUkt2MyW3lxQzZTk0OZmux4b9M7Z6dQ12LmgwxtuMjegmrGxweh62HajjtSowMRAuqaW3n2xrO579zUk9+U7UTYssqVDZ3UJQOkXQwBMarFk6J4wCsbc0iJ8Gf+qCjXL5Lznfam6mSphZ1H28C15+UL5/8JSvbDztc8f31lULFaJZnFtby5JY+fvLeLWY9/y+jffcH+EzVOXSfPliHuOG3Prq1Xcnn3m+xyyhuwShjhYJCs1wmunpLAhqNlnZckuqiwuulkPXJdMWx9ASZcDVFjTh4UdxbcsRp8QxFvXMrVAXtUJnmAUkFyO5/sKaSwprnHWuSO0hNDmJYcxqubcmgyWThUVNcr/ZHbC/A28OwNZ/PRvbO4YEJsl8fZZ8d3+UNAp4NR50PWt6ftulUUxXE7j2udbpbNSnbvA3LmKvDyh5S5Lp0eE+xDca2Hyy3sxlwCyXPg2z9pY6uVIaPJZGFLdgVPf3uU217dRvofv2bJv7/jdx/vZ9OxCibaNqp/m1nq1HWP2zbldZVJDvEzEurn1WMm+WiJFkSnRjle5nT+uBikhLWdrVlKqC+FvM3aeHgHFVW365H83T/BYtLu2HYUlqIFytHj+EPT44wvWuHwayh9x9DfCxgorFbJ8+uPkRYTyILRzmeB7piTwt1v7uDJNUcxWay9umnPbmk3wbFdRIC2u77bT8qjlsDO17UNBSMWeGp5ijKkvLIxh0AfA1eeneD6RaxWrT9y6nng5Vif1Y6ig3w4WlLu+hq6IwSc/2f471z4/hltKp9yxttbUM0NL26lvkUrsxkVHcBFE+OYmhzKlKQwEsO08oIl/97A1pwKINXha+dVNBIZ6I2fsetwJDnCv8eBIlml9egEDI/svrMFUtomSErGRfuSGmjm2O714OUFlcegIgsqjmn/mWwdpIQObvwfjFzY7aXNFiuldc1aj+Tq45DxKpx1E4R3kXjzj4BbP+Xo01dxf+0z8K0fLPiN9n2mDAgqSLb5+mAJWaX1PHndpE7LFnpy3phoksL9eGHDMYC2T9X9rS2TXNdNlnj4fDD4wJEvVZCsKC4orG7ii/3F3DE7BX9vN36sFu6E+mIYfaHLl4gJ8qGsvgWLVaLvjZKv2HRt0+/3z8C0u7Q3euWM1dxq4Wcf7CHA28CT101iclIoIX6dtzacMTyc97Yfx2S2YjQ4dqP6eGXXnS3sUiL82ZxV0e0xWaX1DAvz0/bl1BbCmke1jabW1ragGHlqzbQAVgMU2P5DQMgwLahNnAZhIyBsOKz5P/hwGfxgDUR0/QGgpK4Fq7T1SF7/Vy3YnffL7v8CjP5sOPvf7Prml1y34e9a9vriJ1WgPECoINnm+fXHGBbmx4UOZGc7o9cJbp+VwiOfHCDM39i2Y7W/tQXJDd1kko1+kDJPy2AteVx9cyqKk97ckoeUklvOSXLvQpmfgdDDqMUuXyI62AeLVVJe30J0kGvZ6B4t+A0c+gQ2PqFllpUz1r9WHyGrtJ43bp/G3FHdj1ifMTyM1zbnsregmikO7us5XtnIOSPCTz4gJRTv0zov6bUQJSXcnxU7T9BksuBr1Hd6naOldYyL0MPav8Cm/4C0QPp14BehvacJHWD7td2fj1a28rftrdx9xWKmTDobDN6nXzwqDV5YAO9eB3d+A76hna6hyNYKdriuCHa/q32IDO75zlJ8eBD3mn/AeVPGEbHzGfALh/Me6fE8pfepINnmr1dOpLy+5bSJP864ekoC/1p9hEmJIS5lo3uDr1GPv1HffSYZtLrko19B2WHtB4KiKA5pMll4Z+txzh8XQ0Jo9xmxHmV+pg3v6OJN2BGx7abu9VqQHDkK0m+AbS/CjHshOL53XkfpV9tzK3nxu2xumD6sxwAZYFqKFuxuzal0KEhubrVQXNtMUli7Eond78DKeyEgGiZeA+k3kBwRAkBuRUOnrVnNra2cXfEZv2lcDnnl2vCb8/4Aock9riGx1cJ3O78mvjCIKVM6CZBByy5f+xa8frGWUb7xf20BfHuFtkmXYw8/o92dnfOzHl8f7ANFBDtG3s/5unrY+C8IioNpP3DofKX3qI17NqNjApk10r3bhn5GA+/+YAb/d8k4D63KM7ocKNLeqCXar0e+7P0FKcoZZMWuAmqaWrl9thvDQwDKs6D8sMtdLezsvZKLPDWauivzfqndvt7wt959HaVfNJrM/PzDPSSE+vLrC8b0fAIQ5m9kVHQAW7K7L42wK6hqREoYFt7uzuuedyE4UesRvuU5eO4czvvuapbpv+DEiU7atWWvx/LfeTxu+C+mgHi4/Wu4+jWHAmQAHy89s0ZEsCazBCll1wcmnQMXPQHZa+Hr33R6SFF1E2niOEFZK2HGPRDg2P6mtl7J1c1w4b9g9AXw+S/g4CcOna/0HhUke9jYuCASe6iv6mvh/saeg+TgeIiZqIJkRXGC1Sp5ZWMOE+KDmZLkevYXgMOfab+Odm4UdUf27LHHp+51FJoEU5bBrre0TU7KGeXxLzI5XtnI369KJ8CJOvvpKeHsyKtyaJLdyR7JtkxybSHkbtQ2u13/Djx4GJb8FaNBzyNeb7Lgsznw7vVa8FhyAN65Dt64BGtjNfeZ7ufEFSth2HSnv9aFY6LJr2ziaGn3beY4+2aY8SPY+jzseO20p4tqmvml8X/gHQwz73f49UP8vPA36rU2cHoDXPkyJEyF5Xdq3TWUfqOC5CEgIsCbinoH2ruNWgL5W7UJQYqi9Oi7rHKOlTVw++xk90usMj/TPqiGJLp1mXB/I1560Tu9kjua83PQeWlTO5Uzxsaj5bzxfR63z0phxvDwnk9oZ/rwMBpNFof6JduD5LZpe/uXA1LrKwzaptAZ96C/ZwPXGZ7gu7Br4MRO+OBmeG6m1pHpvD/w2tkfsMp6DiOcaP/W3rlpWsZ3zSEH2tct+qM2EOSzB7WAvh3v4p2cKzJg1v1OlUwJIUgI9Ts5UMTop03mCxmm1UGXZjp8LcWzVJA8BDhUbgEweol2+/To6t5flKKcAV7dlENkoDcXTohz70L1pZC/ze1SCwCdThAV6ENJb5dbAARGa7eV932oZfaUQa+2uZVf/m8PwyP9+cX5o50+f1qKVou8NafnZMvxykb8jXrC/W3dMvZ9CHFnd9oyzRo5hqcNt8BPD8CNy7Vg9f6dMPunZJa3Eh/i63JnmZhgH8bHBzk2otqe6Q1NgfdvhqrctqeWlL1ErS4Epv/Q6TUkhPqeOnXPLwxuWq7VNr91JdSccPqaivtUkDwERPgbqWw0Ye7p9lfsWdpmiSNf9M3CFGUQyyqtZ93hMm6ekeRwu6suHf4CkJDmXqmFXXSQt0uZZJPZSpGzI61n/hi8A+Fb1eXiTPDopwcprm3mn1ena+3UnBQV6MOISH+H6pKPVzYyLNxfuwtTdgSK9pzMIneQEu5PbkWDFqSmngezHoAAbTNhVmk9Ix2ctNeVc9Oi2Xm8quvptO35hmiZXmnVyj9a6iDnO85q3c2G6JvA2/m1JIb5UVDZeGpddGiStkmwuQbevgqaqp2+ruIeFSQPARGB3kgJlY09fPPrdJC6GLLWqOl7itKD1zbnYDTouGH6MPcvdvhz7dZq9Hj3r4V96p7zQfIb3+cy7+/rKHYmC+0XpgXKhz+Dgh0AbDhSxhOrjzj9+kr/+uZgCR/uKOCH80dw1jDXa+ynDw8nI7eqx8RMXkXDyR7J+z7U2rKNv6LTY5Mj/CmvN1Hb3HrK4xarJKu0nlQ3g+SFaVFYJaw/4uDEwPAR2gbBssOw/AdY1/yRIhlG3vDrXXr9hFBf6lrM1DaZT30idiJc9xaUH4X3bgRzL03TVDqlguQhwN4r2aG65NFLoaUWjn/fy6tSlMGrpqmV5TtOcNmkuLbvL5e11MOxtdoAEQ+1jowJ8qW4prn73fqd2Hm8CpPZynvbjzv3gjPu0frRfvtHpJT8cdVBnlmbhcXq3Osr/aeqwcRDK/aRFhPIjxc6PjGvM9NTwqhvMXOwqLbLY6xWSX5VE8PC/bTeyPs+1EaxB8Z0enxKhLa5r+PkvRNVTbSYrW5nkifEBxMZ6M03jtQl241YAEv/Cke+QFewjafMlxMV6togMXuHi/yqTsZvD58Plz8PeRvho7u1yZxKn1BB8hDQNlDEkbrk4fNB7626XChKN74/VkFTq4Wrp7i3yQ6AY2vA0gJprk/Z6ygm2JtGk4W6FnPPB7eTWaSN4X1323GHuhO08Q6EOQ9C9joObPqUrNJ6zFZJhSM/c5QB4Xcr91PdaOKf16TjbXC+zKI9+2a/rdld1yUX1zZjMlsZFuYHJ3ZAVU6XpRZwMkjO6RAkHy3V/s2mRrsXJOt0gnNHR7HhcJlz//an3gkz76c2cjIfWOYRF+LaIDF7j/VT6pLbm3AVLP6TNkXw85+rQLmPqCB5CAgP0DZFOBQkG/21T/OHv7CN8lQUpaOM3Eq8DTrPjJ/P/FzbCT/sHPevZdPWBs6Jsokmk4WcigYmJgRTUtvCNwcd2MTU3pTbISge3+/+Amg/O/qkw4bitlV7C1m1t4gHFqYyLs79f9PRQT4kh/uxNafruuRTOlvs+1BLzoy5uMvj7R0wcstPzbRm2dq2jYx0rbNFewvHRFHXYma7A5sO2wgBi//EN+e8gRkDscGuDfBp65XcWSbZbub9Wi12xsuw4k5VetEHVJA8BLRlknuaumc3eon2qb78aC+uSlEGr+25laQnhridccPSqt21GbWk0wleroqxT91zIkg9XFKHlHD33BHEh/jy5pY8517Uy4eqqT9hRMshfpyQpb1+X3TYUNxS32Lmdx/vJz0hmB/OP72rhKtmDA9nW05llyU3+ZW2IDnEG/av0Ka++nQdoPt46YkP8SWn/NRexkdL64kK9CbYz8vtNc9OjcBo0LEm04mSC5tC20jq2GDXMsnBvl4Eehu6ziTbnfd/WmeP/cu1zXzNXZe0KO5zKEgWQiwRQhwWQmQJIR7q5PlfCCF22/7bL4SwCCHCHDlX6X1BPgaMeh3lDQ5+6mybvqe6XChKR40mM/sLa5ma7ObwEIDs9dBc7fYAkY7sU/ecCVIzbfWj4+ODuGH6MDYfqyDLdivbUS/WziTXGs0Pre8hsKpM8iBwtKSOqsZW7l0wEoPec3mz6cPDqG02k1nceRCXV9mAQSeIr9oGDaXdllrYJUf4kVNxaqb1qAc6W9j5GQ3MHBHOmkM9TN/rRGFNM6F+XvgaXfvgLIQgPtS3+0yydqCWTb78v9qgkdcugDon7/ooDuvxO0IIoQeeAZYCY4HrhRBj2x8jpfy7lHKSlHIS8DCwXkpZ6ci5Su8TQhARYHQ8kxycANET4LCqS1aUjnYfr8ZilUxJDnP/Ypv/AwExWlcZD3Jl6l5mcR3+Rj2JoX5cOzURL73grS2Ob+BrbrXwzo5Cvom5E9/KQ1xi2KoyyYOA/f+R/Xa/p0xP0eqSt3RRl5xX0Uh8qC/6A8u1CXUOfA8kh/uTU1bfFsBKKTnmgc4W7S1MiyK3opHsDrXPPSmqbnI5i2yXGOZHdpmDr5t+ndaGriIbXl6kjbVXPM6Rj43TgCwpZbaU0gS8B1zazfHXA++6eK7SS8IDHBwoYjd6CeRvUdP3FKWD7blVCAFnu9EiC9Amh+Wsh3PuBS/X6hi74uOlJ8TPy6lM7sGiWkbHBKLTCSICvLlgQizLdxTQaHJs898newqpbmxl7OLbIGocD3otp7jayZ7LSp8rsgXJ7gZ4HcWF+JIY5svWLvolH69sZESIHg59CmMvduh7ICXCn9pmM1WNWhu44tpm6lvMjIx2vx7Z7twx0QCODRZpp6immbgQ976PpyaHkl3e4PiHy5HnwW2fgqkBXlnc1oJR8RxHguR4IL/dnwtsj51GCOEHLAGWO3uu0rsiAozOBcmjlmqN0rO+6b1FKcoglJFXyejoQIJ93ayB3PRvLYM2eZlH1tVRTJAPxTWOfc9LKcksqiUtNqjtsZtnJFHXYmbl7kKHzn99cy6jogM4Z0QkzPoxw2QhIeUZLq9f6RvFtc0YDTpCPVDT29H0lHC25VZi7aQuOa+ikUWGXWCqc6jUAk7vcHG0RKtP9mQmOT7El7SYQMdGVLdT6IFM8pxUbTjKd0fLHD8pfjLc8TUYA+D1i9TEXA9zJEjurHFnV8U6FwObpJT29KPD5woh7hJCZAghMsrKnPgHojgkIsDbsT7JdnFngX+UbRKYoigAZouVnXlVTHW31KI8Cw5+AlPvAJ+gno93gTZQxLFMbmFNM7XNZsbEnMzITU4KJS0mkDe/z+uxPnPn8SoOFNZyyznJ2vS0MRfTpPNjarX6+THQFdc0Exvso/1/87DpKWFUN7ZypENte01jKzVNrcxoWKuVGyXPceh6yR16JR+1d7bwYJAMcN6YaDLyqqhpbO35YKChxUxts5lYNzPJaTGBRAR4893RcudODB8Bd6yG8JHwzrWw+x231qGc5EiQXAC0bwaaAHSVWriOk6UWTp0rpXxBSjlFSjklMjLSgWUpzogI9KaiocXxzQg6HYyyTd+zOPaDQlHOdJnFdTSYLExNcTNI3vwk6I0w44eeWVgnnMkk2zftjWmXSRZCcPM5SRwsqmVXfnW357+2OY9AHwOXn2W7UWj051DoQua1bkS2OLf5T+lbxTXNbd1QPK2rfsnHKxsJop6kyo0w/krQObbZLTHUD71OtGWSs0rrCPXzItzf6NF1nzsmCotVss7B6Xv2Ue5xbmaShRDMSY1gU1Z5p9n3bgVGw7LPIWUOfPxD2Py0W2tRNI4EyduBVCFEihDCiBYIf9LxICFEMDAPWOnsuUrvC/c30mqR1DQ5EfCOWgotNWr6nqLYbLP1T3Wrs0VtEex5D866CQKiPLSy00UH+VDR0OLQYITMYi2QHR1zam3nZZPiCfA28Nb3XbeDK61t5ot9RVw9ORF/75Nt7PKTLsNPtNC8d2WX5yr9r6i2qa0biqclhvkRH+J7Wr/kvMoGluq3o7O2akMyHGQ06EgI9SWnwh4k15MaFejxLPikhBDC/Y0Ol1wUVtvrut3/e5yTGkFFg6nbaYVd8g6EGz6EsZfB17+BPe+7vZ6hrscgWUppBu4DvgIOAR9IKQ8IIe4RQtzT7tDLga+llA09nevJL0BxTGSgfeqeEyUX9ul7qsuFogBaPXJ8iK97tYdbngWrWRsM0Itign2QEkrres4mHyyqJTHMl0CfU+tS/b0NXHl2PKv2FlHZ0PnPjre3HsdsldxyTtIpj4th55BrjUbuetv1L0LpVVarpKSmpdeCZNBKLrZmV55yFzOvopFLdZuwho7QSvuckBzuT255A1JKjpTUM9LNSXud0ekEC9KiWHe4FLMDHzLbMskuTttrb/bICAA2ZjlZcmFnMMIVL2glLCt/pLWZVFzmUFNEKeXnUspRUsoRUso/2x57Xkr5fLtjXpNSXufIuUrfc2o0tZ13gHbrRo2oVhSklGzPrXIvi9xUBRmvwLjLISzFc4vrRNtAEQd2ymcW1ZIW03lt9E0zkjBZrHyQkX/acyazlXe2HWf+6Mi2etG21w/25X+WufgVboYqJweTKH2istGEyWIltpfKLUDrl1zRYGqbjAdQXZzHDP0hdOnXaH1/nZAS4U9OeQNl9S3UNLUyMtLzQTLAeWOiqG02k5FX1eOxOeWNCHGy9aI7ooJ8SIsJdG7zXkcGb7j2La1G+f2boETlJl2lJu4NEU6Npm5v1BKoPKam7ylD3vHKRsrqWtzrj7z9ZTDVw6yfeGxdXXG0V3Jzq4Wc8oZTNu21lxodyIzhYby9Ne+06WlfHiimrK6FW2cmn3ZebLAPKyxzkAjY8+5pzyv9z/4BKsbD7d/aa+uX3G7U87CiL9AhYbzjpRZ2KRH+NJosfH9MK+FI7YVMMsDs1EiMeh3fdjN9r7C6iQfe28Xz648xMSEEo8EzIdWc1Ai251TRZLK4fhHfELjxQzD6w9tXQ80Jj6xtqFFB8hBxcjS1C0EyqGyyMuRtz9UySi53tmhtgq3Pa71NYyd6cGWdc3Tq3pGSOqzy1E17Hd08I5n8yiY2HDk1u/X65lySwv2Yl3r6ZuuoIG8KiSA/ZKq2297a821rpW+d7JHce5nkpHA/ooO8T+mXPLVuDcd9RkPESKevZ79jsfqg1sc4NcpzPZLbC/A2MH14GN900i+50WTmX6uPcO4/1/Hl/mLuP3ck79w53WOvPTs1EpPFyrZcN+cUhCTCDR9Acw28c40aYe0CFSQPEaF+RnQCKrqoK+xSSCJEj1d1ycqQl5FbSZCPwfWerLvfhoYymP1Tzy6sC6F+XhgNuh4HimQWaZv20roJkhePiyYy0Js3t5wsm9h/ooYdeVXcPCMJne70W+beBj1h/ka+D1oC1XlwfLOLX4nSW4pttbS9GSQLIZgxPJytOVpdsqn4EGkym+xY10axD7cFyesOlxHobSA6yNuTyz3FwrQosssa2rppWK2SFTsLWPCPdfxnzVEWjY1hzYPzeHDx6FM2rbprWnIYRoOO7454oB1u7ES45g0oy4QPbgazkzHAEKeC5CFCrxOE+Ts5dc9u1BKtw0VTz7VZinKm2p5byZTksE4Dwh5ZzLDpPxA/BZJmeX5xnRBC2NrAdR8kHyyqxddLT1KYX5fHeOl1XD81kbWHS8mvbAS0LLKvl56rpyR2eV5MkA9rmQbGQNW7dQAqqmnGoBOEB/ReoAlayUVZXQs55Q3UZ7yHVQoaU10bvhsX4otRr6O+xcyIqIBe6e9st7Dd9L0deVVc/txmfvbBHmKCfPjfPefw1PVnkRDa9feNq3yNeqYmh7q+ea+jkQvh4v9A9jr49MfgaCtYRQXJQ0lEgJGyOhc+RY5eCtKi9UxWlCGoor6FY2UNrpdaHPxYy6bO/qnTG5XcERPk03MmufjkOOruXD99GDoheGfbcaoaTKzcU8jlZ8d3O3kwJtiH43XA+MvhwMfQUt/lsUrfK65pJjrIB70rH/ycMH249n2z5VgFPpkr2GwdS3R8Ug9ndU6vEySGaTXUnpy015nEMD9GRQfw5JqjXPncZoprmvjXNel8dO8s9/YmOGBOaiSZxXWUOjFavltn3Qjzf63tD1j7F89ccwhQQfIQEhHgYiY57izwCYHstR5fk6IMBvYd7i51tpASNv4bIkbBaNduMbsqOtin2417Ukoyi+sYE9tzXWdssC/njYni/e35vPF9HiazlVvPSe72HG3qXzNMuglaG+CQapM/kBTVNPdq+ze74RH+JPmbidj+N/zqj7PSOovEbu5c9CQlQguOe2vTXnuXpMdhMlv58cJU1v58PlecneDa3SQnzUnVWsE5PX2vO/N+CWfdDBv+Bjte99x1z2AqSB5CIgKMVDS4ECTr9JAyV+u3qG7TKENQRm4lRoOOCQnBzp+ctQZK9mkdLXR9+yM3Jsib4prmLidtFtc2U93Y2u2mvfZunpFMZYOJ/3x7lBnDw04bPnL66/tQ2WCiOWYyhI0A1TN5QCmu7YMg2dSA2PRvPpP3sbjiLQ6ELGCNfjaRbpR4pERoAXZvbdpr7975I9nzyGJ+tmgUfkbP1R33ZExMEOH+Rs+VXIB2F+uiJ2DEQlj1U3V32AEqSB5CIgK8KXel3AK0wSI1+VCZ7dE1KcpgsD23ivSEYLwNjo3PPcXGJyAoHiZc7fmF9SA6yIcWs7XLSZttm/a66JHc0cwR4QyP8MdildzWSdu3juy9mkvrTDDpBsjbCJU5ji1e6VVSSoprmnuvR3JrM2x5Dp5Mh2/+QE3YRC5q+RM/sf6UqLBQt2qJx8QGodcJ0hy4A+IunU7g4+XC970HXnd2agTfHS3v8kOuS/RecM3rWp92VXbRIxUkDyHhAd40tVpoaDE7f/Lw+dqv2es8uSRFGfCaTBb2n6hxrQYxf7sWGJ7zI20SVh9rawPXRcnFoWKtJZSjwYZOJ3jgvFTmpEZwnm1Tk8Ovn34dILSR3Eq/q20y09Rq8Xwm2dKqDcz5z1nw5UMQmQa3f0Xj1e+xXw7naGk9w9wotQC4dFI8q386173Jl4PA7JERlNe3tI2N9xjvQO1D+4kd0FDR8/FDmAqSh5AIVweKAIQNh+BEFSQrQ87u/GrMVulaPfKmf2v1/Gff6ullOaSnqXuHiuqID/ElyKfrzXcdXTopnjfvmI5B3/PbxylBcnCC9mF7j+qZPBAU1drbv3ko0JRS+wD09BTtVn5wAtzyCdy2CobNYGRUAOH+2ntQUrh7QbJeJxjeS5P2BpI5tv7jbk3f60rqIkDCsW89f+0ziAqSh5CIQPtoahdKLoSA4fMgZwNY3ZgCpCiDTIatof/kYU5mkvO3QeYqmH63NuK9H/Q0UCSzqNahTXvuv74WkDHpRqg+Dnmbeu01FccUtU3b81AmeftL8NHd4B0EN3wId3ytvWfYCCGYlqJ9D7mbSR4qYoJ9GBUd4NnNe3axZ4FfBBz92vPXPoOoIHkIifC3B8kuZJIBhi+A5moo2uO5RSnKALctt5LR0YEE+zmebaUqD967AUKTYcYPe21tPYkK7LrcornVQnZ5g8Ob9lwR6G3Az6inuMb2M2fMRVoQpXom97tiTwbJ5Vnw9e+0aZJ3rYdRizttdTjdHiSH+7v/mkPE7JGRbMuppLnVw8kpnU7rn3xsjbqz0w0VJA8hEYFulFuA1uECIGe9h1akKAOb2WJlZ14VU1OcKLVoroF3rtUmW93wIfi6UKbhIUaDjogAY6dt4LJK67FYpcOb9lwhhLC1gbNlkr18YfwVWt/oFg/XWSpOKappRgiICnRzkIjFDB/fAwZvuOTpbju4XDIpnttmJjOtl3sMn0nmjIqgxWxlu7sjqjuTuhgaK6Bwl+evfYZQQfIQEm7PJLva4SIgCqLGqbpkZcjILK6jwWRxfIiIxQwfLoOKo3DtGxA5qncX6IDoLqbuHSpybtOeq06b+jfpRmhthIMre/V1le4V1zQRGeCNlwO15d3a/CQUbIcL/wlBsd0eGuZv5A+XjMPX2PfdIgar6SlhGPU6NvZGycWIc0HoVMlFN1SQPIQYDTqCfb1c65VsN3w+5H0PrU0eW5eiDFT2emSHO1t8+ZB2+/LCf57sCNPPtKl7p3/PHyqqw8dLR3Iv3/qOCfKhpP3rJ0yF8JGq5KKfFdU0E+tuqUXxPlj7GIy9DMZf6ZF1KafyMxqYnBTKht4Ikv3CIH4yZK32/LXPECpIHmLCA4yul1uA9sZvaYH8rR5bk6IMVNvzqogL9iE+xIEOAFv/C9tfhJn3w+Tben1tjupq6l5mcS2jowN7fSRxjO31rVZbr1chbD2TN6m+6/2o2N1pe+YW+OgerZzown/16bj1oWbOqAgOFdVSVufGe3dXRi6CEzuhoReC8DOACpKHGLcGigAkzQSdQZu+pyhnMCklGbmVjmWRj3ylZZHTLoLz/q/3F+eEtql37Tb+SCk5VFTbq5v22l4/2AezVVLe/g5W+vXabV7VM7nfFNc0u9f+bd3jULIfLnkK/MM9tzDlNHNGaq3gNnly+p6dvRWcmr7XKRUkDzGRAd7uZZK9A7TbpaouWTnDFVQ1UVLb0nN/5OL98L/bIWYCXPGCNsZ9AGmbeteu5KG0roWqxlbSehgr7cnXP6UuOShO65az7QUti6X0qbrmVupazK5nkvO3aT3Az7oJRi/x6NqU042LCyLUz4sNvdEvOXYS+EeqkosuqCB5iIlwt9wCtJKLwl3QVOWRNSnKQLQtx4F65LpirZOFdxBc/x4YB15rq86m7p3ctNc3mWTopFfz0r+BMRBeuxCOqI1DfclefuNSTbKpQSuzCEqA8x/z8MqUzuh0glkjI9jo6RHV2sVhxELI+kbNQOiECpKHmPAAb2qbzbSY3fhmGD4fkJDznaeW5XFWq+REtdpceCZZsbOA/Sdq+uz1MvIqCfQxMDq6i2yrqRHevR6aKuGG97Ts6ADUeZCstV8b04vt39peP+j01wcgYiTc+Q1EpMK718GO13p9LYrG3rfa/v/GKd/8ASqPwWXPgE/v//tRNHNTIymta+FISb3nL566SEt6qbs6p1FB8hATEaC1gatwZeqeXfxkMAYM6JKL1YdKmPe3tSpQPkMcLq7jZx/s4YpnN/PO1uOez6Z0YntuFVOSQtF1trHNatWmixXugitfhtj0Xl+Pq6JtgVBJu0xuZnEtccE+zg1IcVF4gDcGneh86l9gNNz2GYxYAJ8+AN/+WRtvrPSqohoXR1Jnr9NKZKb/8GTffKVPzE6NAHppRLW9FZwquTiNCpKHmIgAbaCIW0Gy3guSZg3ooSJ5FQ2YrZJDhbX9vRTFA97ddhyjXsf04WH8+qN9/PzDvTSZeu/WYGWDiazS+q5LLdY9Boc+gcWPQtoFvbYOTwjyMeDrpT8lk5tZVNcnm/YA9DpBVKB3p1P/APAO1EpVzroJNvwNVv4ILK19srahyv6BJSrIiUEizTXw8Y8gPBXOe6SXVqZ0JS7ElxGR/r0zotovDOKnqH7JnXAoSBZCLBFCHBZCZAkhHurimPlCiN1CiANCiPXtHs8VQuyzPZfhqYUrrokIdHM0td3w+VCRBdX57i+qF9g/BGSV9cKtKaVPNbdaWLGzgPPHx/D6smn85LxUVuwq4IrnNpNb3tArr7kjT6u373SIyIGPtGBu0k1wzn298vqedHLqnRYYtZgtHCur7/UhIu3FBHc+0KSN3kub1jb/Ydj9NrxzjZrI14uKapsJ9zfi4+XEJtMvHoK6Irj8v9rkRKXPzUmNZGtOhedHVIM2fa9wF9T3QqZ6EOsxSBZC6IFngKXAWOB6IcTYDseEAM8Cl0gpxwFXd7jMAinlJCnlFI+sWnFZhG3qXpkngmQYsNnkcnuQXKqC5MHu831F1DabuX5aIjqd4CfnjeLV26ZSVNPExU9vZPXBEo+/ZkZuJUa9jokJwac+UbQXPr4XEqbBRYOnN2x0kHdbuUVWaT3mXh5H3VH7IL1LQsD8h7RgOXs9vLpU2xipeJzTPZIProQ978Ccn0HC5N5bmNKtuaMiaG61sjOvFzbNp56n/XpMtYJrz5FM8jQgS0qZLaU0Ae8Bl3Y45gZghZTyOICUstSzy1Q8JSJQK7dwO5McNQb8owZsXbJ9qqAKkge/97blkxzuxznDT/ZinT86ik/vm01yuD8/eCODv32Zidli9dhrbs+tZEJC8KmZtvoyeO8GbXjCtW+BwYlb1f0sJsiHIluQnGnftNdH5RZwcjS2Q7XkZ98MN3wAFdnw0nlQmtn7CxxinJq2V30cPrlf24sy71e9uzClW9NTwvHSi96ZvheTrrWCUyUXp3AkSI4H2t9TL7A91t4oIFQIsU4IsUMIcUu75yTwte3xu9xbruIuP6MBP6PevZpk0LI+w+dpGZ8BuNHG/vUdK6vvk01eSu/IKq1jW24l100bhuiQtU0M8+PDe87h+mnDeHbdMW55ZZv7H/7Qyjv2nahhSvv+yGYTfHAzNJTBdW9rG84GkZhgX0rrtKl3h4pq8TboSA7367PXjw32odFkoa7F7NgJqefBss/BYtI+mKjvYY8qrmlq29DZLYsZlt+p/f1f+bJWFqP0G39vbUT11weKT06w9BSdDkaeB8e+Va3g2nEkSO7sfmLH/zsGYDJwIXA+8DshxCjbc7OklGejlWv8SAjR6ZZYIcRdQogMIURGWZmqielNbo+mths+HxpKofSQ+9fysIr6FnQC6prNvTPKU+kT727Lx0svuGpyQqfP+3jpeeyKCfz9qonsyKviov9s5EiJe7WsK3efoNUimTdKm3KFlPD5z+H493DpMxB3llvX7w8xQd60WiSVjSYyi+sYHROIQd93+7Y767DRo7hJsOiPWrux49/3zsKGoOZWC1WNrY5lktc9Bvlb4aInICyl9xen9Oj6acPILm9g/ZFeiJPaWsHt8Py1BylHfkoWAInt/pwAFHZyzJdSygYpZTmwAUgHkFIW2n4tBT5CK984jZTyBSnlFCnllMjISOe+CsUpEe5O3bNLmaf9OsBKLqSUlDeYGBun3U5WJReDU3OrheU7C1g8NqatdWFXrp6SyIp7Z9JitvDoqoMuv6bFKvnv+mzGxQWdLO/Y/hLsfB1m/wwmXOXytftT+4Eeh4pq+2TSXnv2VmNFzgTJAGMu1tpN7nm3F1Y1NNk3UMb01P4tez1890+t68gg/Xd/JrpgQiyxwT68tDHb8xcfvkBrBadKLto4EiRvB1KFEClCCCNwHfBJh2NWAnOEEAYhhB8wHTgkhPAXQgQCCCH8gcXAfs8tX3FFRIA35XVullsAhCRC2IgBFyTXt5gxma1MS9aCHNXhYnD66kAx1Y2tXD9tmEPHj4sL5u55I/juaDk7j7u2sWX1wWKyyxu4Z94Irbwjez188SsYtRTO/Z1L1xwI7Jnc/SdqqGgw9emmPehmoEhPjP4w9lI48DG0qp7nnmD/oNJtJrmhHFbcpQ16Wfq3PlqZ4ggvvY5bZyazKauCg55uceoXBglT4ajql2zXY5AspTQD9wFfAYeAD6SUB4QQ9wgh7rEdcwj4EtgLbANeklLuB6KBjUKIPbbHP5NSftk7X4riqIgA77aNbW4bPh/yNg2ovqb2euSxcUEEeBtUJnmQenfbcRLDfJk5Irzng21unpFEqJ8XT6056vTrSSl5bn02w8L8WDo+Bipz4MNbtUDhihe0mr1Byp5JXndYu0Xbl5v24GQ/3m7bwHUl/TpoqYXMzzy8qqGpuFb7sNFldwurFT7+oXbb/apXBuSo9aHu+qnD8DPqeXljjucvnroIinZDveq/AA72SZZSfi6lHCWlHCGl/LPtseellM+3O+bvUsqxUsrxUsp/2x7LllKm2/4bZz9X6V8RAUYqG0xYPFH4P3w+mOodqmGSUvLVgWLqHd284yL7B4CIACMjogJUkDwIZZfVsyW7kuumDut84l0X/L0N3DlnOGsPl7G3oNqp19ySXcme/Gp+MHc4BnPDyQ1j17876MfvRgZ4oxOwMUvbFT+mD3skg1Y7HuZvdD6TDJA0G4KHqZILD7FnkrscSb31Oe12++I/QcyEPlyZ4qhgPy+umZLIJ3tOUOrK91R3Ri7Sfs36xrPXHaQGb2pEcVlEgDdWqU0Vc1vKHEA4VHKxNaeSu9/cwZXPbuZ4RaP7r90Fe4/kiABvRkaqIHkwem97Pgad4OopnW/Y684t5yQR5GPgqW+znDrv+fXHiAgwcvVZMbDibig7DNe8DmHDnV7DQGPQ64gI8Ka+xUxssA8hfsY+X4O9DZzTdDpIv1bbdV9b5PmFDTElNc0E+Rjw9zac/mThLlj9CIy+EKb9oO8Xpzhs2axkzFbJG9/nefbCMRO19q6q5AJQQfKQZN8E5ZHNe76h2i50B4Jke7BaUNXIJc9sZHNWL/R65GS5RXiAkZFRAZTWtVDbPHDKQZTutZgt/G9HAeeNiSYq0ImBBzaBPl7cPjuF1QdLHK7ZO1hYy/ojZdw9NQyf96+Fw5/BksdODs05A9hvr/f1pj272J6m7nVn4nUgrbDvA88uagjSeiR3smmvpQ7+dzsERMGlTw+aQTlDVVK4P4vHRvPW1jyaTB5s2abTaSUXx77VWgAOcSpIHoIiArQsktu9ku2Gz4eC7T2Okc0tb8DboGPVj+cQGeDNza9s4/XNuR7vY1xhC/7D/LUgGeCYyiYPGqsPllDZYOK6aYk9H9yFZTNTCPQ28PRax2qT/7vhGOOMJdyeeSfkbtSmvk2/2+XXH4jst9fT+rge2S46yIcSV28NR4zUphzuflf1THZTcW0X0/Y++zlU5cIVL2obuJQB7845w6lubGX5zgLPXnjkedBcDScyPHvdQUgFyUNQuCczyaAFyVYz5G3u9rDcigZSIvxJifBnxb0zWTA6kkc+OcDDK/ZhMntuWlpFg4lAHwPeBn1bkKxKLgaPd7cdJz7ElzmprreCDPbz4rZZyXy+r5jDxd1/eMuvbKRm35cs9/od+pZauPVTberbGcYeGPX1pr221w/yoaLBRIvZxaxX+nVQdgiK9nh2YYNUQVUjD6/YS2axcx0OOp22t+c92PsezP0lJM/y4CqV3jQlKZSJCcG8sjHHs8NFRthbwamSCxUkD0GRng6SE2eAwUdrl9WN7PIGksO1ndKBPl68cPMU7lswkve253PDi1s8NvSjvL6lraQkMdQXo16n2sANEnkVDWzKquC6qYnondiw15nbZ6Xgb9Tz9NpuapOlZP/yx3jZ8Ff0YUlw11pIOset1x2o7G3gxvRjuQVAaa2L3+fjrwC995DfwCel5IPt+Sz593e8uy2ft7ccd/hck9lKeX3LqZnkumL47EEYNhPm/qIXVqz0FiEEd8xOIbu8gXVHPNiNwjdUu3OTpYJkFSQPQUG+Box6HWWeCpK9fCBxerd1yWaLlfzKRlIiT7YT0ukEPz9/NE9dfxb7C2u49OmN7D9R4/ZyKhtMhPtrJSUGvY7kCD9VbjFIvLc9H71OcPUU10st7EL9jdx8TjKr9hZ2fifB3ELz8ntZeuI/HAyajdedX0OIYz2ZB6NL0uN4YGEqIyID+uX1o22BmdMDRex8Q2H0Utj34YBqOdmXSuuaufP1DH65fC/j4oJITwhme26lU+dL2aFH8o7XtQ5FlzwF+k428ykDWttwke883A4udZF216auxLPXHWRUkDwECSEIDzB6riYZtJKL0gNdfkOdqG6i1SJJCT+95+bF6XH8756ZAFz1/GY+3dNxoKNzKupNhAec3L0/UrWBGxRMZisfZuRzblpU1z1cnXTnnBR8DHqe7ZhNri+DNy7FZ/87PGm+HN+b3gbv/gke+0pimB8/XTTKqZZ6nmQPzFxqA2eXfj00VgzJ28Cr9hay+IkNbMwq53cXjeXdH8xg4ZhoDpfUUdPk2IeG06btWcyw4zUYca5W960MOl56HbfNTGbzsQoOFLqfZGqTqlrBgQqSh6zwAKPnyi0ARi7Ufv3+qU6fzilvACA5ovPG9OPjg1l532zGxwVz/7u7OFrSfR1pdyoaWtrqrgFGRgZwvLKR5lYP7gBWPG7NoRLK601c78aGvY4iAry5acYwPt59glzbv0GK98GLC5CFu/glP2H/qPsYGR3ssddUOmcv9yiucWNy3siF4B8Je97x0Kr6x18+P8RDy/eyfEcB+ZWN3W5erm40cf+7u7jvnV0khfnx2Y/ncMfsFHQ6wZTkUKSEnXmOTZg8bdrekS+grhCm3OH216T0n+um9cJwkZiJEBgLBz7y3DUHIRUkD1ERAd6eDZJj07UftJufgn3/O+1pe5Cc0kWQDBAZ6M0fLx0PuL7RzmKVVDaYiPA/mUkeERWAVWobB5WB651tx4kN9mHeqCiPXvcHc4fjpdfx7LosbXPpy+eD1cwnZ7/MB83TuGfeCI++ntK5IB8DfkY9xTVu/NzRe8GEa+Dwl9DoeJnBQNLQYuaFDdl8kJHPgx/uYc7f1jLz8W/58bu7eHNLHkdK6to2Ya09XMriJzbwxb4iHlw0iuU/nNm2GRngrMRQDDrhcMnFyUyyLUje/jIExcOoJZ79IpU+FeyrDRf5dE+h6x1kOhICJi/T6pJLMz1zzUFIBclDVESAN+V1Hiy3AFjyuLb5Y+V9p+1Azy1vIMDb0NZ+risxbt6SrW40YZWcmklWHS4GvPzKRjZmlXOtBzbsdRQV6MMN04eRt2sN1reuhKA4Wu9Yw1/3+DItOYzJSaEefT2lc0IIYtxpA2eXfh1YW2H/cs8srI/ZP6w/ed1ZfPmTOfzx0nFMTgpla04Fv/t4P4uf2MDZf1rNVc9tZtmr2wnx8+LjH83i/oWpGPSnvmX7GvWMj3e8Lrmophk/o55AbwNUHIPstTD5NlWLfAa4fVaKbbhIrucuOvVOMPh2eYd4KFBB8hAVHmCkoqHFsz2KDUZtQplfGLx3IzScHBaSU9FISoQ/oocG9aF+Xhj1Okpc3AFf0XBykIjdiMgAhIBjpSqTPFC9vz0fAVzjgQ17nblvRAUvG/5KhQiH21bxyTFJYU0z98wf/NP0BpPoIB+K3Cm3AIidCNHjB22Xi9xybdroiMgA0mKCuOWcZJ6+4Wy2PLyQ9b+Yz9+vmsjisdHUt5i5Z94IPr1/NuPjuy4Hmpocyp78GofKyYprm4gJ9tF+Dme8AjoDnH2Lx742pf8MC/fj/LExvL31OI0mDw0B8Q+Hs26EvR9oXVCGIBUkD1GRAd60WiS1TR6eqBMQBde+BQ1l8OFtbbvQc8rru6xHbk8IQVSQt8vZJnsJSbj/yUyyj5eehFBf1QZugDKZrXyQkc/80VHEhXQyCcxdx7cS/vH1NPtEclnDQ5ywBPPfDccYHR3IgtGeLe1Quhcb7OPyB+BTpF8HJ3ZA2RH3r9XHcsq1n0PJEX6nPC6EICncn6unJPK3q9L58idzeWhpGt4GfbfXm5ochslidagzUFuP5NYm2PUWjLkYAmNc/2KUAeWOOSm24SInPHfRGfdq7+Nb/+u5aw4iKkgeoux9hD3WBq69+LPh4ich9zv4+re0mC2cqGoiJdyv53PBrVuy7UdStzcyUnW4GKg+31dEaV0LN5+T5PmL52+Ht66EgChab/6EUsL4wesZHCmp5+55w3u8s6F4VnSw9r3t9uCDCdeA0A/KbHJ2eQOxwT74GT1T4mAvF9rmQMlFSU0zMUG+sH+FNlFNbdg7o0xJCiXd08NFwkdoH6YyXoaWofceqoLkIcoeJFf0RpAMWqZnxo9g6/NUb3oVq+SUHsndiQ7ycbkmuaItk3xqkDwiMoDssnosnpxKpLhNSsnLG3MYEenPPDcm7HWqIAPeugL8I+DWVcQmDOeqyYkcLKolPsSXi9PjPPt6So9ig30wWyXlDW7+3AmM1jpd7H0frJ6b1tkXcsobut3A7KzwAG9GRPqTkdt9hwuLVVJS16Jlkre/BBGjIXm2x9ah9D8hBHfMGU5OeQPfZnpwuMjMH0NzDex603PX7MBqlR6dvOspKkgeouyZ1nJP9kruaNEfIWUekesfYpLIapu215PoIB+Xp3JVNJjQCQjx65BJjgqgxWzlRJWb9ZCKR2XkVbHvRA3LZqV4tn9vwQ5483KtPv62VRAcD8C980fgZ9Rz74IReOnVj7++Zm8DV+JOhwu79Oug9gTkbnD/Wn3I00EywLSUMDJyK7vNHpbXt2CxSsZyDAp3wtQ7tA4Gyhll6fgY4oJ9+PVH+/hyf5Fn9h0lToVh58D3z2q9tT2kxWxh3eFSfvPRPmY8toZ3tuZ57Nqeot4lhqgIT4+m7ozeAFe/RoMxkueNTzDcx7FbNdFB3tS3mKlvcf6bsbzeRJi/8bQOCW0dLspc77+seN7L3+UQ7OvFlWcneO6iJ3a2C5A/g+CT104M8yPjt+dxw7Qzd7LeQBYTZJ+654EPq6MvBO9g2D14Si6qGkxUN7Z6PEiekhRGbbOZI6Vd/3yz90ieVLwcvPy0DxnKGcdLr+OFW6YQHuDNPW/t5M7XM8ivbHT/wjPvh5rjcPBjty5T09TKyt0n+NE7O5n86Dfc9up2Pt51ginJoaRGB7q/Tg9TQfIQFeZvRCd6OUgG8AvjlcQ/EywaCf70DjD3/HptbeBcGF9bUd9yyqY9O9UGbuDJr2zk64PF3DB9GL7G7jcnOaxwF7x5GfiGwK2rTgmQ7fyMBlWL3E/sQyw80svVywfGXw6HPhk0tZLZtn7xwx0sPXPU1OQwALZ3U3JRXNNEEPXEHF8FE64GHzVA50w1Pj6YT++bxW8vHMP32RUsfmIDz68/RqvFjXKGUUshfCRs/g84mZ0urWvm9c253PTSViY/upoH3tvN1uxKLk6P49XbprLjd4t49sbJzBoZ4fr6eokKkocovU4Q5m/s3XILmy0NcTwT9FPI3wpf/LLH46MCtTfSUhfeSCsaTKdt2gOt/CIiwKiC5AHktc256ITgFk9t2KvO1zLIPsFaiUVI77STU1wXHuCNXifasppuS78eWhu1QHkQaJs86mDpmaMSw3yJDvJme07Xm/eKapq5Sv8dOkuz1v9WOaMZ9DrunDOc1T+bx+zUCB7/IpOL/rORHXkuDuHR6eAc2wyE3O8cPs1ilVz69CYe+eQARTVN/GDucFbcO5Ntv17IY1dMYEFaFD5eHkqS9AIVJA9h4f4enrrXhZzyBgoTlsLsn8GO1zqdyNeeOwNFKupPHUnd3gjV4WLAqGtu5f3t+VwwIZbYYA+0fbNa4KN7tFZFN38MIaqcYiDS6wTRgd4ub8w9TeJ0CBsOe97zzPV6WW55A3qdIDHMsU4/jhJCMCVZq0vuSnF1EzcZvkEmTNV6TStDQnyILy/eMoUXbp5MXXMrVz73PQ+v2Et1owsJsvTrwC9Cm6zroKzSeopqmvnTZeNZ8+B8frUkjbOHhXp2D0ovUkHyEBYRaOz1ILnRZKa4tpmUcH8497cQPxm++FW3I2Wjg7Qg15V+qhX1ptM6W9iNjArgWFmDUxsZCqoaedOTE4wUAD7MKKC+xczts1M8c8HN/4G8jbD0r1rLImXAig72camUqlNCwLgrIHcjNFR45pq9KKe8gWFhfr2yaXRqUiiFNc2cqO683juwaDPDRRFCZZGHpMXjYlj9s3n8YE4KH2QUsPCf69lwpMy5i3j5wvS74ejXUHrIoVP25FcDMHNEuJMrHhhUkDyERQT0fibZPl0qJdIfdHq4+D9af86vftPlOX5GA4E+BqfrFlvMFupazF2Ovh4ZFUBNU6tTJSb//uYov1t5wDMbjRRAu/322uZcJieFMikxxP0LFu6Gb/8MYy6BSTe6fz2lV8W40eKxU2MuBmmBI1947pq9JLsXOlvYTU2x1SV3UXIxtWwFtSIIxl7WK6+vDHz+3gZ+c+FYPr1vNkG+Xjy66qDzF5lyhzaqevPTDh2+u6CaIB+Dx0uM+ooKkoewiADvtuEbvSW3okMNXsx4refinncge12X50W7MFCksm0kddflFuD45r3mVgtf7tdGcWYWq64YnvLNoRKOVzZy+ywPZJFNjbD8Tq0X8sVPqpZWg0CMLZPskdZUALHpEDwMDn3qmev1EqtVktuLQXJaTBAB3ga2d1ZyUVvI5Jbv2R56gbbhURnSxsYFcfusZI6W1nOkxMn3Nv9wOOsmrUe5A6Oq9+RXk54YMmjKKzpSQfIQFh5gpNFk8dyc907YN6qc8sYw75daHeGnP9HGo3bClWxT27S9bsotAIfHU3+bWdrWhi6zSAXJnvLKxhziQ3w5f1y0+xf7+rdQcRQue05r+aYMeDFBPjSatLs+HiEEjLkIjq2FloH7fVpS10xTq6XXgmS9TnB2UminQ0VkxmvopORo4tW98trK4LNkfCw6Aav2FDp/8jn3andvtj7f7WHNrRYyi+tITwhxbZEDgENBshBiiRDisBAiSwjxUBfHzBdC7BZCHBBCrHfmXKV/tPVKruu9bHJOeQNRgd74e7cbwerlq2X9qnJg3eOdnhcV5O30QBF76UhXmeTYYB/8jXqOOZhJXrn7BJGB3sQE+ZBZXOvUWpTO7T9Rw9acSm6dmYTB3brMw19qo1LPuQ9GLPDMApVeZ9+YW+KpumTQSi4sLXB0teeu6WE5Zbb2b70UJANMSw7lcEndqZuyLK1Yd77OButEfKJUvb6iiQz0ZsbwcFbtc2HgSNhw26jqV7r9YHqgsAaLVZLuibK6ftLju5QQQg88AywFxgLXCyHGdjgmBHgWuERKOQ642tFzlf4TaQsmy3qxLrnL24spc2HSTdou2aK9pz0dYyu3cGb+vD2T3FVNshCCEVGOdbioaWplbWYZF0+MY2xckMoke8grm3LwM+q5dqqb3SfqS2HljyB6PCz8vWcWp/SJkwNFPBgkJ04H/8gBXXJh75Gc4uEeye1NsfVL3pHXLpt8+HP09cW8aTmPGE90klHOGBdNjCO7rIFDrry/tY2qfqvLQ3YdrwYgPXHw9uR2JJUzDciSUmZLKU3Ae8ClHY65AVghpTwOIKUsdeJcpZ/Y+wlX9GKQ3O0I1sWParfIP/2x1sKrneggH8xWSaUTbWoqGrrPJAOMdLAN3Ff7izFZrFw6KY60mECOldUPyLnyg0lpbTOf7ink6skJBPt6uX4hKWHlfVoG48qXwND1/29l4HGnxWOXdHoYfYG2677Vg9f1oJzyBny8dEQH9l5N8KTEELz0gm3t65K3v0yTXxxrrWe1DXNRFIAl42PQ6wSf7XOh5CJhSo+jqvcU1BAf4ts2+2AwciRIjgfy2/25wPZYe6OAUCHEOiHEDiHELU6cq/STk6Ope6fcoqaplYoGE8ldBcl+YbDkcW1KWofapugg56fuVdSb8Dbo8O9metuIqACKa5upa26FmhPa5kFL62nHrdxzguRwPyYmBJMWG4TZKjnmYC2zq8wWK00mS88HDlJvbcnDbJXc5u6GvYyX4ehXsOiPEDXGM4tT+owr39sOGXMJmOohZ33Px/aD3PIGksP9e3UDk4+XngnxwSfrkgsyIGc9mfFXYkWngmTlFGH+RmaOCGfVXhdKLkDLJtcch23/7XQKn7Zpb/BmkcGxILmz7+iOfxsGYDJwIXA+8DshxCgHz9VeRIi7hBAZQoiMsjIne/cpLrFnknurDVxuZ5v2Ohp/JaQuhm//BFV5bQ/beyWX1jn+RlpebyIiwLvbkcNpIVau0a+F1y+GJ8bBG5fCk+mw6UloqtZes7aZzccquGRSPEII0mK0efK9XZf8+08OsPTJDU6VmAwWza0W3tp6nIVpUe5tXCo7Al/9FkYs1Pp1KoOOj5eeUD8vz2aSQSvh8g4asNP3csobPD6OujNTk8PYW1BNc4sJPnsQAmNZH3o5Bp3o9i6bMjRdNDGWvIpGDhS68P42aolW6vTVr+GF+XD0m7ZgubLBxPHKxkG9aQ8cC5ILgPbzXROAjrn5AuBLKWWDlLIc2ACkO3guAFLKF6SUU6SUUyIjIx1dv+IGb4OeIB9D7wXJFQ4EyULAhf8CBHz2s7ZvsLZbsjWOr62ioaXTkdSYW7RaxfdvYsGq2fzN60VE7QmY9yu46lVt+MTq38O/xsLnv2Tdlm1ICZdOimtbv1Gv69W65LK6Fv6XUUBuRSMZeafvTh/sVu4+QWWDyb3hIWYTrLhT2/h52bOq3dsgFhPs69mNewAGo/amnfl5l7d/+0urxcrxysZe62zR3pTkMFotkqK1/4Wi3bD4Txyv1xMd5IN+kLbhUnrP+eNiMOgEn+51oeRCp4PbPodLn4WmSnj7SnhlCeR8x56CaoBBvWkPHAuStwOpQogUIYQRuA7o+FF9JTBHCGEQQvgB04FDDp6r9KPe7JWcU96AEDCspxGsIYmw8HeQ9U3byOrIAG+EcK5usaLeRJi9/ZuUkLsJPrkf/pEK798Ex7cgJ9/GFa2P8sz4D2DBwzD+Crj1U7j7Oxh7CWS8wlWbL+HdwKcY0bgXpMRLr2NkVECv9kp+d9txTBYrRr2OVa78sBrApJS8vDGHtJhAzhnuxtSltX+Goj1wyVMQGOO5BSp9LibI27Mb9+zGXKS9WR//3vPXdkNBVRNmqyQlIqDXX2tKUiih1BKb8TdImg3jr6S4prkt8aAo7YX4GZmdGsFnrpZc6A1w1o1w3w4t4VWdB69fxIjPb2Cy7ggT4s/wcgsppRm4D/gKLfD9QEp5QAhxjxDiHtsxh4Avgb3ANuAlKeX+rs7tnS9FcUVEgHevdbfIKW8gLtgXH6+ua4TbTLtLG1n95UPQWIlBryMiwJtSp4LkFsL9bbcTv38GXrsA9i3Xsks3LoefZaK/4G/UhqeTZWvH1CZ2Ilz+PHm3bOUZ86WcJQ/Cq0u1W0gHPyEtNrDXyi1MZitvbslj3qhIzhsbxef7ijBbzpxNgpuyKjhSUs8ds1O6LYXpVu4mrSTm7Fu0QEgZ1GKCfZ0eFuSQkeeBwWfAdbnIKdf2M/RFJjnU38ifAlfgZa6HC/8BQqggWenWRRPjKKhqYk9BjesXMRhh6h3w411w/mOE1B1lufEP+P/vBi25MUgZej4EpJSfA593eOz5Dn/+O/B3R85VBo6IQGOvZUidmi5lH1n9wjxtZPXlzzk1UERKSXmDSWv/VpoJa/6o7Xa/8iUwnrqGkZEBXU4Z+ijLwpOWa7j6nieIyfkItjwHH9zM+RP/zYraKCob2mWrPeTzfUWU1bWw7KpkGk0WPt9XzNacSmaNjPDo6/SXVzblEBFg5OL0ONcu0FwLH98Doclw/mMeXZvSP2KCfKhoMNFituBtcOBDtKOM/lqgnLkKlv51wJTkZPdBj+Q2BTtY2rqaN+WF3BSRhk5KimqaOTctqvdfWxmUFo2Nxksv+GxvIZPcLY/w8kXO+CFLVyfwm+jvuCD/A/jvXDAG2r4fhW23mtD+LHQnfz/rJzDzPne/HI9yKEhWzlwRAd6U15V7/LpSSrLLG7hskhPNTOwjqzf+C9IuJDoomhPVjgXJ9S1mTGYrEX46+Phu8A7Qgm7j6W9KI6MCWH2oBJPZitFw8maKlJJPdhcyIyWcmIgwiLhDG7/50kIWHP4/IvkTmcW1zBzhueBVSsmrm3IYHuHP3NRIWsxW/Ix6Vu0tPCOC5Ir6FtYdLuWeeSMcu6PQma8ehpoCWPal9v9VGfRigm0bc2tbSOypHMtZaRdpQXLhTu3u1ACQU95AiJ8XoR7+gH0aqwU+f5AW7wj+XnMZU4vriA/xpanVojLJSpeCfb2YmxrJZ3uLeHjpGLc7sBRUNXGiUU/V2fdB+m9hx2vaCGspAXnqr9J68vcRozzw1XiWCpKHuHB/b2qbzacFjO6qbDBR12zuuv1bV+b9SmvL9tHdjE95np21jv1gt9dVzyx6U2spd/VrEND5BtCRUQFYrJLcigZGRQe2Pb7/RC3Z5Q3cNXf4yYMN3nDly3j9dx7/9Hqew0XneDRI3nm8mj0FNfzx0nHodAJfo55FY6P5Yn8xf7x0PF7uTqXrZ18fLMEqtdt5Lsn8TGtWP/tnMGy6Zxen9Bv7UIvi2mbPB8mjzgedQSu5GEBBcnJ4H2SRd74BhbtoPP8Z6lf6kZFXiU6nDRiJVYNElG5clB7LmsxSduVXMTkpzK1r7c6vBtA6W/gEw6wH3F9gPxnc78CK2yICbQNFGjxbl2zvbOH07UUvH7j2LfDy5dbjv6a1oZoWc8+9gysaWhgj8hhz9HkYdwWMu7zLY0dGadnIjkNFVu4+gZdesHR87KknRI6G8//MXP0+Iva/4tzXA9on5C9+Bc/M0EpB2nltcy6BPgauPDuh7bGLJsZR3djKpizPZ/j72hf7i0kK92NMbGDPB3dUXwaf/BhiJsD8hz2/OKXf9MrUPTu/MEieowXJrmxE6gU55Q29X2rRWAlr/g+SZhM2/QZignzYllPZ9nesMslKd84bE43RoGPV3iK3r7Unvxpvg47RMS783B9gVJA8xLUNFKnzbIcLew2e05lkgOB4uOYNgpsLecLrGUprmno8pbK2nn96PY/FOwQu/Ge3x9p7lR5rFyRbrJJP9xYyf3QUwX6nT4MTU24nw+cclhQ/3+kY7S5JCV/8UhuWUp0HLy+CY98C2jCFL/YVce2URPy9T97UmTsqgkAfA5/ucf+HVX+qbjSxOaucpeNjnd+wJ6U2ibGlDi5/QdsUopwx7AGbx9vA2Y25GCqyoOxw71zfCU0mC0U1zb2/aW/N/2n1+xf8HaHTMTUljO25lW1DW9QgEaU7gT5ezB8Vyef7itzu1b+noJrx8cGD/k4oqCB5yIvopYEiuRUNGHSChFAXb/ElzSTr7N9wnn4X+g2P93h47J6nGavLo27RP7RMUjf8jAbiQ3zJajdBb2tOBSW1LW29kU8jBGtH/54qGYhcfieYGnv+GqTUNiFuewHOuQ9+tA2CE+GtqyDjFd7ckotFSm6dmXzKad4GPYvHxvD1wWKHsugD1eqDJZitkgsmuNCubddbcPhzWPh7iB7r+cUp/SrIx4Cvl753MskAaRcCYkB0uWjrF9+bg0RO7IAdr8P0e9q+X6Ymh1JS28L23Ep0AiID1SARpXsXToylpLbFrV79ZouVfSdqBv0QETsVJA9x9kxyWZ1ng+Sc8gYSw/zc+iRpnXIHH5jnEbfnqe7f7E7sZGzWSyy3zMF/4sUOXXtEVMAp5RYrdxXib9Rz3pjoLs9JSkzkZ633IMoPw9e/7f4FpIRvHoEtz2hvXIv/pPWDvv1LGHEurPopMd8/yqK0iE5rMi9Kj6Wu2cyGI4O35OKL/cXEh/g63yezKldrBZg8B2bc2ytrU/qXEIKYYJ/eaQMHWh/txGkDYvpejiOTR91htcJnP4eAKJj/UNvDU2x1pV8fKCEiwPuMyOopveu8MdF4G9zr1X+kpJ7mVuugH0dtp75rhrjYYF/C/I2uTdvpRk55I8nh7m3IiQ7y5XfmZZQFjYOP7jmtnheA1mb4+IfUe4XxT/3tDreTGhkZwLGyeqxWSYvZwuf7izh/fEy3HRjGxASx0TqB7NRlkPGyNtmrM1JqY7Y3PQlT7oAlj59sReUTBNe/x9Gk67mZVTzW+jcwNZx2idkjIwjx8xq0g0Vqm1vZeLScpeNjnCu1sFrgox9qbYEue1ab6KSckWKCfChyoJTKZWMuhuK9p4y77w/2ILnXNu7tekPr5LH4T9rPF5vRMYEE+hiobzGrUgvFIf7eBs5Ni+LzfcVYXCy5sE/ac7uV3ACh3oGGOKNBx91zh/Pd0XIycis9ck0ppa1HsnvtukL8vJAGH94f/pg2ivi9G6Cp+tSD1v0FyjJ5K/rneAc4viN3ZFQAza1WTlQ3se5wGXXNZi7toV1danQAOgGfRtypbSb75D6trU1H6/8K3/1DG3xxwT9O69UqdXrur7mBZ33vIuzEt9oYz9pTg2EvvY6l42NYfbCEJtPgK7n49lApJouVpRNiez64ve+fhuObYenfIGRY7yxOGRCSwv3Yf6KWf3x1mNrmVs+/QJpt6EzmKs9f2wnZZQ3EBPmcsu/AYxor4Zs/QNIsmHD1KU/pdYLJSaGA2rSnOO6iiXGU17ewNafCpfN3H68mxM+r50m7g4QKkhVuPieJiAAjT3xzxCPXK6ltoanVQkqEe98kQgiig7w52hwE17yhbXxbcZd2exEgfxtsfgrOvpWNchLhTvQgbetwUVbPJ7sLCfc3MmtE9yOTfbz0JEf4c6CkGa58RatL/uiek+sB2PAPWPcYTLoRLnqy00zoluxKMovrCD/3fsT170NlNrx47mlTiS6aGEejycLaw6UOf10Dxef7iogO8uYsZ7IJxfu1DPyYiyH9ul5bmzIw/GzRKM4fH8PTa7OY97e1vLwxx7M1+GEpED2h3+uSc8rre6/UYs0fbZv1Tv8wDjA1WbV/U5yzIC0SXy89n7nY5WJPQTXpCSGuT1cdYFSQrOBnNHDPvBFsyqpgW4772eSTNXjuD36ICbLVLSbN1MoWjn6lZY/tAWpQApz/ZyrqTYQHOB8k78mv5ptDJVw0MRaDAzV7Y2KCOFxSB5GjYMlfIHstbHlWe3LTk/DtozDxWrjkqS5LBV7dlEOon5eWuR61GG7/CoReyyjvfBOqj4OUTE8JIyLAOOhKLhpazKw/UsbS8bGON6U3t2gfgHxCtA8XZ8gPWKVrUUE+PHX9WXx632zGxQXz6KqDLPznej7aVeD27vo2Yy6G41ugvv8+aOaUN7jW5acnO16HHa+eslmvI3uQrDLJiqP8jAYWjoniy/3FmC3Wnk9op6HFzJGSOtLPkFILUEGyYnPj9CQiArx5YrX72eS2Gjw3M8mgvZGW1No2FU69EybdBBv+Dm9fDZXH4NKnwTuQioYWwgMc370d5m8kzN/I65tzaTFbucTByYBpMYHkVTTS0GKGyctg9IVa66UvfgWrfw/jr4RLn9XGbHciv7KR1YdKuH7asJP1zzHj4QdrIDJNK+H49wT4axKGNy7iqdD3CTn8Po15O7T660Fg7eFSWsxWlo53oqvF2j9D6QHt/6d/9xl95cwyISGYt+6czpt3TCPY14ufvr+HC5/ayPojZUh3+xyPuQiQ2lCaflDdaKKqsdXzPZJ3vqm1SExdDOc90uVhkxJDuH5aYrcbkhWlo4smxlLRYGJLtnNJs/0narBKmHSGbNoDFSQrNr5GPT+cP4Lvsyv4/phrtUh2uRUNGA064jxwi8+eSZZSatnFC/+pTdHK2whTfwDD52GxSiobTEQ4OfJ1ZGQAVY2tJIT6cvawEIfOSYvVNsYcLqnT1nPJU+AbpvVBHnOJ1tNX33Xt4eubc9EJwc3nJJ36RGAM3PE13LEaLvyXFmxbTEyrWsVfdP/F79Vz4S9x2kCSdY8PmCEJnfliXzERAd5MSXawRjzrG9j0Hzj7Vm1amjIkzUmN5NP7ZvPkdZOob2nl1le2ceNLW8ksrnX9olFjIWx4v5Vc9Epni11vwyf3w8jz4Jo3tamgXTAadDx2xcS2O2eK4oj5o6PwN+qdvotp37R3prR/AxUkK+3cOH0YUYHePPHNEbcyONllDSSH+7k9/x20ILnRZKGuxaw94OUD170D5/4OFv0foGVrrBKnMskAI6K0N65LJ8U5XD+VZpsglFlUpz3gHw7Xv6ON077qlW4D5IYWM+9n5LN0fEznNYJ6L61t1dQ74KIn4M5vEA8VcLXX0zwX9XuY/VPwj9Bqntf13Du6PzSZLHybWcr546LRO/L/v+aEVmYRNUYrp1GGNJ1OcOmkeNb8bD5/uHgsmcV13Pl6hus/j4TQSi5yNpy+6bcPtAXJnuqRvPtdWPkjGD4frn1b+3moKB7m46XnvLHRfHmgmFYnSi725NeQGObr9HvxQKaCZKWNj5eee+ePYFtOpVvZ5NyKBo+1O4oK0r7ZStv3Uw2Mgbk/B6P2GhUN2rTAMCczyWNsWeGeulq0lxDqS4C3gcPts1vxk2HBr7UgtxsrdhZQ12xm2awUh19PZzAwYeJknjgxlpqZD8Gtn2qbAtc/DjvfcPg6fWX9kTKaWi1c4EhXC0sr/O92rR75mjfAeGbshlbcZzTouG1WCj9fPJqCqqbTRsg7ZcwlYG2Fo197boEOyilvQK8TJIZ64N/2nvfh4x/C8Hlw/bsqQFZ61UUT46hubGVTluO9+nfnV59RWWRQQbLSwXXThhET5ONyNtlilRyvaPRY5iQmSHsjKK7petiJfVqgMxv3AK6ZksjKH81iVLTj8+WFEIyOCeRQcZ1Tr2W1Sl7dnEt6QrDDpR12F6XHYrJYWX2wRMuMXfwkjFgIn/4EjvT9G393vthfRKifF9NTHCi1WPN/kL8FLvkPRKT2/uKUQWfuqAhA+/DlsrizITAW9q/w0Kocl13eQGKoL0aDm2+1ez+Ej++BlDlw3btaS0xF6UVzR0UQ6G3grS15DsUCZXUtnKhuOmP6I9upIFk5hY+XnnsXjGB7bhUbnfgEaVdY3YTJYiXFQ5nkaFuQ3N1krop6LZMc4eQtHh8vvUu7cNNiAsksqnXqQ8SGo2VklzVw26xkp1vjnJUYQnyIL5/usdWH6b3gmtchehx8eCuc2OnU9XpLc6uFNYdKWTw2pudOIZmfae37pt6p1V8rSicSQv0YEenPhqNuTJ7U6bS7L0e+hLLDnlucA3LKGtyvR96/HD66S+uFfP176o6L0ie8DXp+dO5IvjlUyid7eq5N3muvR1ZBsnKmu3ZqIrHBPjyx2vlssqc3qtiD5OJug2RbJtnJcgtXpcUGUdtspqjG8W4Tb36fR0SANxdOiHP69YQQXJQey6ascqpspSV4B8KNH4JfBLxzDVTmOH1dT9t4tJz6FjNLJ/TQ1aIqV7ttHDsJzv9LXyxNGcTmjopka3YFza1u9FCeca+Wfd34RJeH7Mir5FCRG5sEO5BSkuPuUKUDH8HyH8Cwc+CG99tKzBSlL/xgznDOGhbC71ceOLXksRN78qvR6wTj4oK6PW6wUUGychpvg54fLRjJzuPVTmdwPB0k+xr1BPkYuv0GrWwwoRMQ4tdHQbJt895hB0suCqoa+fZwKddNTXT5tuvFE+MwWyVfHmg34S8wBm5aDlYzvHUlNLjXlcRdX+wvJsjHwMwREV0fZG6BD28DCVz9Wrc78xUFtCC5xWxlqzs93P3DYfJtsPcD7UNaJ37xv7385fNDrr9GB24PVTq4Ev53h7aZ94YPVICs9Dm9TvCPq9NpbrXw64/2dZs0211Qw6joQPyMvTBZsh+pIFnp1DVTEokP8XU6m5xT3oC/UU9koOeCn5hgn24zyeUNJsL8jY51U/CA0bYg+ZCDrane25aPAK6f7vqY5XFxQSSH+53ekidylHYLtvYEvHutNmSlOy11sPsdeOMyePl8bWLXsW/B1ODy2gBMZiurDxazaGxM9x8EvvoNFO6Cy57VJqIpSg9mpIRjNOjY4E5dMsDM+7X+5ZuePO0pq1VSUNnE8coevn+ckF2ubTZ0KZNcX6ZlkOMna3eMvFULN6V/jIgM4Bfnj+abQ6V8tOtEp8dIKdmTX31G9Ue2U0Gy0imjQcd9545kd3416w47/uZkny7lyZGU0e0HinSior6FcP++y0gG+XgRH+J7sg1cN1otVt7bns+C0VHEh7i+2UYIwUUT4/j+WAVldR3+LobNgCtfgoIMWPEDsHa4LW21QNYa7U3376laqUNVrpaB3vhvePNyeHwYvLQIvvk/7dgW57oJbD5WTm2zufsBIvuXw/YX4Zz7bEMeFKVnvkY901PC3A+Sg+Jg0g2w6y2oPXXkbkldMyaLlRNVTU5PGetKbrkWcLu0iXnn62BpgUuf0UqrFKUfLZuVwpSkUP7wyQGKOykzzKtopKap9YzrbAEqSFa6cdXkBBJCfZ3qdJFb4fkRrNH20dRdcHYktSeMiQ10aMjB1wdKKK9v4aYZST0e25OL0+OwSq2DxOkLuhiW/g0yV8EXv9SGjZQcgK9/C/8aC29doY30nnQ93P41/HiXNuXvoeNaycbM+wGpZdneugL+mqQFzRmvOjS45It9xQR4G5id2kWpRXkWfPJjSJgG5/3Brb8HZeiZmxrJ0dJ6Cqub3LvQrJ9oHxq/f/qUh49XaAGt2Sqd2mvQnZzyerwNOmKDnGzVZjFr33cp87Q7RYrSz/Q6wd+vTsdksfLwir2nxQN7ztBNe6CCZKUbXnod9587kr0FNXybWdrj8SazlYKqJo+PYI0O8qa0rgWLtfNgraLB1OfNy0fHBJJd1kCLufvNRG9vzSM+xJe5oyI98pqpUQF8vOtE5x9apt8Fsx6A7S/BfybBczNhy3PaLdtr3oCfH9WGlAybzt+/PsyP391FYZNem9x13h/gzm9sQfMKmPljMDfDqp/AO9dqt3+7YLZY+fpgMQvHRJ0ctd1eaxN8cAvojXD1qz32k1aUjuaN1r5/3M4mh6XAhKsg45VTavjbl1nke6jkQtu05+/8UKWjX0FtAUz7gUfWoSiekBLhz6+WpLH2cBkf7ig45bnd+dX4eulJPQMnOzoUJAshlgghDgshsoQQD3Xy/HwhRI0QYrftv9+3ey5XCLHP9niGJxev9L4rzk5gWJifQ9nk/KpGLFbpsUEidjFBPliskoqGzksuyutb+qyzhV1aTBBmq+RYade1vMfK6tl8rIIbpg/zWL30DdOHsfN4Nd91taFy4R9g8jLwj4Klf4cHj2gTAcde2rZJ7mBhLc+uO8YnewpZ9K/1vLwx5+QtZu8AGLkQznsE7t6gZaez18Fz58DR1Z2+5NacSqoaWzsvtagvhY/uhtIDcMWLEJzg/l+CMuSkRgUQE+TDhqNuBskAs38GrY2w9bm2h/KrTmaoPVWXnF3uYvu3bS9CUDyMWuqRdSiKp9x6TjLTUsJ49NODp9zV2ZNfzYT44J5bfw5CPX5FQgg98AywFBgLXC+EGNvJod9JKSfZ/vtjh+cW2B6f4v6Slb5kzybvP1HLY19kdtuGKdfTI1htouy9kjsZKNJitlDXbCaiH8otgG5LLt7ZehyDTnDNlESPve4N04cRH+LL377KxNpZZl2ng4v/DXeu1jLL/uGnHfL3rzIJ8vFi1f2zmZIcxqOrDnLZs5va+ly2EQKm3w13rdOC7revgs9/qWWG2/l8XxG+XnrmjYo6+WBtIXzxEPx7Ihz6FBb+HlLPc/vrV4YmIQRzR0Xw3dFy92uGo9K08qStL0BzDaBlj2ODfTDohEeCZLPFqg1VcjZILs+C7LXaB91uRtwrSn/Q6QT/uCodi5T8arlWdmEyW9lfWEv6GbhpDxzLJE8DsqSU2VJKE/AecGnvLksZSC4/K55rpiTwwoZsLnjyO7Zmd95qrK39Wy9kkqHzgSKVtr7BfV1ukRzuj9GgI7OLNnDNrRb+t6OA88fHeLTTh7dBz88WjWL/iVo+76w2uQdbsytYe7iMH84fwfj4YF5bNpWnbziLktoWLntmE3/45AB1za2nnhQ9Fn7wLUz/IWz7L7x4rlbvjDZh8asDxfx/e3ceJXdd5nv8/VT1vq/pdHoJCdlZQocQ9iCgCAiCoA4KHmdRRifuXBWduc6ox+ty76gwoAziNqIyDgoiIjsDcVRISFg7IQlZupN00mvS+1rf+8evqlOpVKeruqu70tWf1zl9UvXr3+9X3/qeTufJt57v81yybA7ZGX441AAPfxpuW+lt0jv1eli3AS68JRFvX2axtUvK6eofHs1/nJQLb4GBw15qEt7q8fzSHKqLsxMSJO/t6GM44OIPkjfcA750OPODkx6DyFSoLc3hC1cuZ/32Vn75QiNvHOhicDiQkvnIEFuQXAU0hj3fGzwW6Vwze9nM/mBmp4Qdd8DjZvaimd08ibFKkqT5fXzr3Sv52d+tYSgQ4K/u/gtffOBVOiOCqV2tPRRmp1Oc4NSH4zUUCXXbm+50izS/j8Vz8sYMkn//ShOH+4a4cRJl38ZybV0VSyvy+dfHtzEUx6qac45vPrqVioJM/vq8k4AjVTOeuuUibjx7Pj/9827e9u3nePS1pqPTa9Kz4IpvwI2/hp5WuPti+Mv32birldbuQa5fMAAProPb62DTz7wOZx/fBNfeCWWLEjwDMhtdsKgMn8Gz2ybRfS9kXp2Xi//n78FgL43tvdSW5FBTkpOQnOQJ1Ysf7PHKM664BvLmjH++SJLcuKaW8xeV8rXf1/Pwq15Z0lSsbAGxBcnRkikjP+fdBMx3zq0E/g14MOx75zvnVuGla6wzs7VRX8TsZjPbaGYbW1oSkHcmCXfh4nIe+9RaPnzhAu57oYG3fftZHgtrbrFrojl44yjLy8BnRG0o0hrqtjfNK8ng5SVvHaND173P72FheS7nLjw23WGy/D7js29fyq7WHn61sXH8C4Ke3NLMpoZDfOqtS47ZYFeQlc5Xrz2V33z0PIpy0vnIvZv48H9s5JmtzTy7reXIl1vJn97+O1rnng+P3kr+r9/H7Rnf4+InroTX7vfaTH/yZS/lo3jyFT1EQopyMlhZU8Szk928F3LhLdDbytCGH9PcNUBNcQ61JTkJWUneOZEg+dX/8la3tWFPTnA+n/HN608H4N+f3UlpbgbVxRMvcXoiiyXpaS8QnlRZDRzV0cA51xn2+BEz+56ZlTnnWp1z+4PHm83sAbz0jeciX8Q5dzdwN8Dq1avj64Us0yYnI41/fMcKrl45j8/d/wp//7MXufK0ufzLO09hd2sPZ09BUJjm91GWl3ncleTpzkkGLy/515v2enWaw4L01/cfZnPDIf73VSsSWi863KXL53Dm/GJue3I719VVe6kOxzEScPzfx7aysCyX95w59ua5utpifvfxC/jx/+ziO09s58ktY1U1+Ttu8s/nn3ruZaHfj527Ds79OORXTOJdiRzf2sXl3P70djp6Bif/idX886D2POxPt5POt6gtzSEjzUdH7xCd/UMUZE28Csuu1m4KstIoiXWMzsEL90DFaVBz9oRfV2S6VBfn8E9XreALv3mVlTVFU/ZvXbLFEiRvABab2QJgH3AD8P7wE8xsLnDQOefMbA3eCnWbmeUCPudcV/DxZUDkpj6ZgU6vLuJ3H7+Au5/byW1PbeeP271mElOxkgxe171oDUVCFS+StZIMXnvq8xYdef2fP99AZpqP61dFy0pKDDPj85cv473//md+8qfdfPQtJx/3/Ac272PbwW6+d+OqcXcgp/t93Lz2ZN5VVz3Oqtr5bB34JIsqCqAw8f85Eol00dJy7/fNjlauXjlv8jdcewtp917Pdf711JSsJSP4d6OxvZdT5k18I9Lu1l4WlOfFHjg0Pg8HX4WrvuttmBWZAW44q4Y9bb2ce3Lq/v4fN0h2zg2b2ceAxwA/8CPn3Otm9pHg9+8C3g181MyGgT7ghmDAXAE8EPxFkQb8wjn36BS9F5lm6X4f6y5exBWnzuXW37zKC7vaR1s2J1pFQVbUXMG27kEy03zkjrOSOhWWjVa46OK8RV4Tje6BYX67eR9Xr5xHUc7Urm6vWVDCJcvm8P3/3sH719RSmBN95WtgeITvPLGN06oKj98RL0J5fmYMmw6L4xixyOSsrC6iMDud57a1JCZIPvlSWvNX8NHDD5Fb+OWEBcm7WntYs6Ak9gte+AFkFsLp753wa4pMNzPj1iuWJXsYUyqmGjPOuUeARyKO3RX2+A7gjijX7QRWTnKMcoJbWJ7HfR8+h/qmTk6ZVzAlr1FRkMnG3e3HHG/tHqQ0NyMpH/WU5WVSlpdxVBm4Bzfvo2dwZEo27EXz2bcv5crb1/P9Z98c85fVvX9pYN+hPr55/ekp+5GYzA5+n3HBojKe296Cc27yP89mPD3nA7y36wu4Pb8nc8m7gMnVSu4fGmHfob7YP1Xrbob633r5/BlT80mciExM6lV+lqTw+YxTqwqnLAibW5BFR+/QMXWa23oGkpJqEbJsbsFohQvnHPf+ZQ8rKgs4Y5rK4SyvLOCalfP48f/s4kCUdrpd/UPc+cwOzl9UOnbLaJEZ5KIl5RzsHOCNg9Ery8TrycBqdvtqsfXfpiDDT1FO+qSC5N1tcW7a2/RTCAx5QbKInFAUJMuMEGoo0hyRl9zeM0hpEjbthSybm88bB7oYCTg2NRxi64Eubjpn/rSu2H7mbUsJOMftT28/5ns/WL+L9p5BPvf21P5ITGaPC5d4/9mbdIvqoIaOfp4ovRFatsAbjwQrXPSNf+EYdrXEESSPDMPGH8PCi1UqUeQEpCBZZoTRhiJdR6+WtnUPUpqbxJXkygIGhgPsaevh58/vIS8zjXeekYBcyTjUlubw/jW1/OeGxtH6rOCVx7tn/U6uPG1uyhZ6l9mnsjCbJRV5CSkF55yjsb2XpporoHgBPPN/mF+cOalayaHybyfFEiRv+wN07tMqssgJSkGyzAijDUXCUgqcc7R2DySl/FvIsuBGxT/vbOPhV5q4tm4eeZnT3072Y5csJjPNx/97/I3RY3c8vYOB4QC3XLZ02scjMpXWLi5nw64OegeHJ3Wf9p5BegZHqC4pgLf+CzS/ztVDj7G3o5eRaG3fY7CrtYc5+Zmx/R544QdQUA1LLp/Qa4nI1FKQLDNCtNbUPYMjDAwHkppusWhOHj6D7z65ncHhADeenZwGGuX5mXzoggX8/pUmXt17mMb2Xn7+/B7eu7qak8vzkjImkaly0dJyBkcCPL/z2M288Wjs8NIqaktyvE53J13IRfvuJmekK2pd9ljE3FSpZRvsehZW/w34p/8/1iIyPgXJMiMUZKeRmeY7KkhuC3XbS2K6RVa6n4XlebR0DXDm/GKWV05NdY9YfGjtQopz0vnWY1v5zhPb8JnxiUsXJ208IlPlrJNKyEr3TTrlIrRBr6Ykx6tPfPk3yBju4tNp99PQFn/KhXOOnS3dLCyPIUje+EPwZ8CqD8b9OiIyPRQky4xgZsc0FGkNdttL5koyHEm5mK6yb2MpyEpn3cWLWL+9ld9s3sdfn3cSlYWp2SpUZresdD9nLyid9Oa9xtEgOfj3ZO6p9Jz6AW7yP8nhPS/Hfb+DnQN09A6NNhoa00A3vPQLWHEt5JXH/ToiMj0UJMuMUZGfddRHoKGV5LIkloADryTV4jl5XHlaZVLHAXDTOfOZV5hFflbauF34RGayi5aUs7O1Z1Kb7BrbeynLyyQn40i6Q9ZlX6KbbJa/9DWvXXQc6psOA7BivHrxr/4KBjphzYfjHrOITB8FyTJjVBRm0RweJPecGCvJ71ldwxOfuYis9Onv+hcpK93PvR86m19++Jwp7/gnkkxrl3grsM9tn/hqckN775FV5KC0/DJ+kvF+5nduhK0Px3W/+v1eY6Flx+s8GhjxNuzNPQ2qz4p7zCIyfRQky4xRkZ/Jgc5+XHB1J7SSXJKrYDDcwvI8Tq2aeEtdkZng5PJcqoqyefaNiQfJjR293qa9CJvmXMce/3x47B9hKPYNfFuauqgtySE/K3qLeACevwua6+GCT3t50CJywlKQLDPG3MIs+ocCdPZ7ZZ9auwfJz0ojMy35K7giMr3MjLVLyvjTm20MjQTivn5oJMD+Q/3UFB8bJFeV5vN190E4tAf+fEfM96xv6mR55XFWkdt3wlNf9Uq+nXJd3GMWkemlIFlmjDkRZeDaegaTno8sIslz0ZJyugeG2dxwKO5rmw71MxJwUVeSa0tyeLR3GcNL3gHrvw2d+8e9X8/AMLvbelhROcanOIEAPPQJ8KfDVd/RKrLIDKAgWWaMyFrJbd0DlCrVQmTWOm9RGX6f8ey25rivbewIK/8WIRQ47z7zCxAYhif+edz7bT3QhXPH2bS36Sewez1c9lUomN6unCIyMQqSZcaoKPBWjUNd99q6B5O+aU9EkqcgK526miL+uL017msbIsu/hQkFyW8Ol8N5H/OqUTQ8f9z7bWnyNu1FTbc4vBce/xIsWKu6yCIziIJkmTEqjkm3GKBU6RYis9pZC0p4fX8n/UMjcV3X2N5Lms+i1hIPBcmN7b1wwWcgvxIe/byXMjGG+qZOCrLSqCqKuJ9z8PCnwY3A1bcrzUJkBlGQLDNGVrqfopx0DnYOEAg42nsGlW4hMsvV1RQxHHC8vv9wXNc1tPdSVZyN33ds0FqYk05BVpq32pyZB2/7CuzfDC/9fMz7bWnqZHllARYZBL/yK9j+OFz6JShZENcYRSS5FCTLjBJqKHKob4iAQ0GyyCx3Rm0RQNyb9xrbo5d/C6kpyRlNyeC090D1Gnjqy9Dfecy5IwHH1qauY/ORu5u9FejqNbDm5rjGJyLJpyBZZpRQQ5FQjWSlW4jMbnPys6guzo4/SO7oi7ppL6Q2PEg2gyu+CT2tcP/feqXcwuxu66FvaITllRFB8iOfhcEeuOYO8KlUpchMoyBZZpRQQ5HW7hOj256IJF9dbTGbGzpiPr97YJj2nsGoNZJDakty2NveRyAQbE1dtQre/jXY/Ue44ywvAO72qmqENu2tCA+S6x+C+gfhos9D+dK435OIJJ+CZJlR5hZm0dI1QHOXt3lPdZJFpK6miP2H+0cr34ynMbhCPF66xeBIgINdYfc8dx18YjPUfQA2/BBuOwOe+To7GptI8xmLK/K883rb4fe3eK2nz//kRN+WiCSZgmSZUeYUZBFw8MaBLkA5ySICdcG85JcaY1tNPl75t5BQAN3Q1nv0Nwoq4ervwrrnYfFb4dlv8LcvXsdnCp8hk2CFjcf+EXrb4Jo7veYhIjIjKUiWGSXUUKS+qROfQVGOgmSR2W7FvAIy/L6Y85JjWUkeDZLbe6OfULYY3vsf8KGn2e6q+Ie+u+HOs+Cpr8DLv4ALPgWVK+N5GyJyglGQLDNKqKHI6/s7KcnNiFq+SURml8w0P6dUFcQVJOdnpVGYPfYq77yibHx2JKAeS1vRqVzf90X+sPIOyMiH9f8KZUth7efieQsicgKKKUg2s8vN7A0z22Fmt0b5/lvM7LCZvRT8+lKs14rEI7SS3NI1QGmu8pFFxFNXU8wr+w4xNDJ2w4+QhvZeaopzjq1pHCYjzUdlYfbYK8lBW5q6AKPgtCvg75+DG34J7/9PSM+K9y2IyAlm3CDZzPzAncAVwArgfWa2Isqp651zZwS/vhLntSIxKc3LHF09VmULEQmpqy2ifyjA1qaucc9t7Og7bqpFyFFl4MZwpB11Afh8sOxKNQ0RSRGxrCSvAXY453Y65waB+4BrYrz/ZK4VOYbfZ5QHK1qoRrKIhIQ2720eZ/NeIOBobO897qa9EC9I7jvuOfVNncwtyKJEm4hFUk4sQXIV0Bj2fG/wWKRzzexlM/uDmZ0S57UiMQvlJauyhYiEVBVlU56fOW5eckv3AAPDgdhWkktzaO0eoHdweMxz6vd3srwyP97hisgMEEuQHC1py0U83wTMd86tBP4NeDCOa70TzW42s41mtrGlpSWGYclsVRHMSy5TuoWIBJkZq2qLxm0q0jha/m38IDl0TuMYq8n9QyO82dJ9bDtqEUkJsQTJe4GasOfVwP7wE5xznc657uDjR4B0MyuL5dqwe9ztnFvtnFtdXl4ex1uQ2SYUJCvdQkTC1dUWs7utl/aewTHPaYgjSB6vDNyO5m6GA+7YdtQikhJiCZI3AIvNbIGZZQA3AA+Fn2Bmcy24TdjM1gTv2xbLtSLxmlsYDJKVbiEiYepqioDjNxVpaO/FzEvPGM94QXJ9tHbUIpIy0sY7wTk3bGYfAx4D/MCPnHOvm9lHgt+/C3g38FEzGwb6gBuccw6Ieu0UvReZJbSSLCLRnFZdiN9nbG44xCXLKqKe09jeR0V+Flnp/nHvV5yTTl5m2pi1kuv3d5Kd7md+ae6kxi0iJ6Zxg2QYTaF4JOLYXWGP7wDuiPVakcm4aEk571tTwynKAxSRMDkZaSybm3/czXuN7b0xbdoDL8+55jhl4LY0dbKsMl9NjURSlDruyYxTnp/J1687PaaVIBGZXepqi3ip8RAjgah7xGns6I0pHzmktiR6QxHnHPVNncpHFklhCpJFRCRl1NUU0z0wzJst3cd8r39ohAOd/THVSA6pLcmhsb2XQETQve9QH139w8pHFklhCpJFRCRljDYViVIKbt+hPpwj5nQL8M4dGA7Q0j1w1PH6/cFNe0r7EklZCpJFRCRlLCjLpTA7PWpecjw1kkOO1Eo+OuWivqkTM1g2V41ERFKVgmQREUkZZkZdbdFxg+R4V5Lh2DJwW5o6WVCaS05GTPvfRWQGUpAsIiIppa6mmG3NXXT1Dx11vLGjj8w0H+VxlI+sKs7G7NggWZv2RFKfgmQREUkpdbVFOAcvNx4+6nhDWy/Vxdn44ijZlpnmp7Ig66ggubN/iMb2PuUji6Q4BckiIpJSVgY770Vu3muIo0ZyuJpghYuQrU1dACyvVD6ySCpTkCwiIimlMDudRXPy2Nx4aPSYc47G9vhqJIfURjQU2TLajrpw0mMVkROXgmQREUk5q2qL2NzQgXNefePDfUN0DQxPaCW5tiSHg50D9A+NAF75t+KcdCoKYs9tFpGZR0GyiIiknLraYjp6h9jT5q0AN7b3AfGVfwupLfWu2dvh3WvLgU5WzCvATO2oRVKZgmQREUk5o01FGr285FC6RE3xxHKSQ/cYHgmw9UAXy+dq055IqlOQLCIiKWfxnHxyM/yj9ZJHg+Q4WlKHjNZKbutlZ2sPg8MBVbYQmQVUBV1ERFKO32esrDnSVKSxo5finHTys9LjvldpbgY5GX4a2vsoDm7aU41kkdSnlWQREUlJdbVFbGnqpG9whMYJln8Dr4tfqMJF/f5OMvw+Ti7PS/BoReREo5VkERFJSXU1xQwHHK/tP0xDey+nVU28ZFtNSQ4Nbb0MDI+wuCKPjDStMYmkOv0tFxGRlHRGcPPext0d7Ovom1Bli5DQSvIWtaMWmTW0kiwiIimpLC+T2pIc/vBaE8MBN+F0C/CC5L6hEfqGRlihIFlkVtBKsoiIpKy62iJe2XsYmFj5t5DwAFsrySKzg4JkERFJWXU1RaOPJ7OSHJ6qoZVkkdlBQbKIiKSsutpiwCsJV1mUNeH7VBd79ZWrirIpzIm/jJyIzDzKSRYRkZS1vLKAzDQf5fmZpPsnvi6Ule6nqiibU9RERGTWUJAsIiIpKyPNx5oFJWSn+yd9r+/duIrinIwEjEpEZoKYgmQzuxy4DfAD9zjnvjHGeWcBfwH+yjl3f/DYbqALGAGGnXOrEzBuERGRmNx105mYTf4+K8Pym0Uk9Y0bJJuZH7gTeBuwF9hgZg855+qjnPdN4LEot7nYOdeagPGKiIjEJTdTH5qKSPxiSdBaA+xwzu10zg0C9wHXRDnv48CvgeYEjk9EREREZNrFEiRXAY1hz/cGj40ysyrgXcBdUa53wONm9qKZ3TzWi5jZzWa20cw2trS0xDAsEREREZGpEUuQHC2Ty0U8/y7weefcSJRzz3fOrQKuANaZ2dpoL+Kcu9s5t9o5t7q8vDyGYYmIiIiITI1YErX2AjVhz6uB/RHnrAbuM29nRBlwpZkNO+cedM7tB3DONZvZA3jpG89NeuQiIiIiIlMklpXkDcBiM1tgZhnADcBD4Sc45xY4505yzp0E3A/8g3PuQTPLNbN8ADPLBS4DXkvoOxARERERSbBxV5Kdc8Nm9jG8qhV+4EfOudfN7CPB70fLQw6pAB4IrjCnAb9wzj06+WGLiIiIiEwdcy4yvTj5Vq9e7TZu3JjsYYiIiIhICjOzF8fq4THxHp0iIiIiIilKQbKIiIiISAQFySIiIiIiEU7InGQzawH2JOGlywC1z556mueppzmeHprnqac5nh6a56mnOZ4e8c7zfOdc1AYdJ2SQnCxmtnGs5G1JHM3z1NMcTw/N89TTHE8PzfPU0xxPj0TOs9ItREREREQiKEgWEREREYmgIPlodyd7ALOE5nnqaY6nh+Z56mmOp4fmeeppjqdHwuZZOckiIiIiIhG0kiwiIiIiEkFBMmBml5vZG2a2w8xuTfZ4UoWZ/cjMms3stbBjJWb2hJltD/5ZnMwxznRmVmNmz5jZFjN73cw+GTyueU4gM8sysxfM7OXgPH85eFzznGBm5jezzWb2cPC55jjBzGy3mb1qZi+Z2cbgMc1zAplZkZndb2Zbg7+fz9UcJ5aZLQ3+DIe+Os3sU4mc51kfJJuZH7gTuAJYAbzPzFYkd1Qp4yfA5RHHbgWecs4tBp4KPpeJGwZucc4tB84B1gV/fjXPiTUAXOKcWwmcAVxuZuegeZ4KnwS2hD3XHE+Ni51zZ4SVytI8J9ZtwKPOuWXASryfac1xAjnn3gj+DJ8BnAn0Ag+QwHme9UEysAbY4Zzb6ZwbBO4DrknymFKCc+45oD3i8DXAT4OPfwpcO51jSjXOuSbn3Kbg4y68X8RVaJ4Tynm6g0/Tg18OzXNCmVk18A7gnrDDmuPpoXlOEDMrANYCPwRwzg065w6hOZ5KlwJvOuf2kMB5VpDsBRSNYc/3Bo/J1KhwzjWBF+ABc5I8npRhZicBdcDzaJ4TLpgG8BLQDDzhnNM8J953gc8BgbBjmuPEc8DjZvaimd0cPKZ5TpyFQAvw42Dq0D1mlovmeCrdAPwy+Dhh86wgGSzKMZX8kBnFzPKAXwOfcs51Jns8qcg5NxL8WK8aWGNmpyZ5SCnFzK4Cmp1zLyZ7LLPA+c65VXhphuvMbG2yB5Ri0oBVwPedc3VAD0qtmDJmlgG8E/ivRN9bQbK3clwT9rwa2J+kscwGB82sEiD4Z3OSxzPjmVk6XoD8c+fcb4KHNc9TJPix6X/j5dtrnhPnfOCdZrYbL+3tEjO7F81xwjnn9gf/bMbL4VyD5jmR9gJ7g582AdyPFzRrjqfGFcAm59zB4POEzbOCZNgALDazBcH/jdwAPJTkMaWyh4APBh9/EPhtEscy45mZ4eW9bXHOfTvsW5rnBDKzcjMrCj7OBt4KbEXznDDOuS8456qdcyfh/R5+2jl3E5rjhDKzXDPLDz0GLgNeQ/OcMM65A0CjmS0NHroUqEdzPFXex5FUC0jgPKuZCGBmV+LlwvmBHznnvpbcEaUGM/sl8BagDDgI/DPwIPAroBZoAN7jnIvc3CcxMrMLgPXAqxzJ4/wiXl6y5jlBzOx0vA0gfrzFhV85575iZqVonhPOzN4C/C/n3FWa48Qys4V4q8fgpQX8wjn3Nc1zYpnZGXgbUDOAncDfEPzdgeY4YcwsB29f2ULn3OHgsYT9LCtIFhERERGJoHQLEREREZEICpJFRERERCIoSBYRERERiaAgWUREREQkgoJkEREREZEICpJFRERERCIoSBYRERERiaAgWUREREQkwv8HxNlSCgV720sAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot between noised and smoothed ndvi\n", "fig, ax = plt.subplots(figsize=(12,5))\n", "ax = df.plot(ax=ax)\n", "ax.set_title('ndvi_2010-2012')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Save result**" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "NDVI_MA.to_csv('output/3/no_noise_data/no_noise_ndvi.csv', index = False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.10" } }, "nbformat": 4, "nbformat_minor": 4 }